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OPTIMIZATION  OF  ACOUSTIC  SOURCE  STRENGTH  IN  THE  PROBLEMS  OF 

ACTIVE  NOISE  CONTROL* 

J.  LONCARIC'i'  AND  S.  V.  TSYNKOV^ 

Abstract.  We  consider  a  problem  of  eliminating  the  unwanted  time-harmonic  noise  on  a  predetermined 
region  of  interest.  The  desired  objective  is  achieved  by  active  means,  i.e.,  by  introducing  additional  sources 
of  sound  called  control  sources,  that  generate  the  appropriate  annihilating  acoustic  signal  (anti-sound). 
A  general  solution  for  the  control  sources  has  been  obtained  previously  in  both  continuous  and  discrete 
formulation  of  the  problem.  In  the  current  paper,  we  focus  on  optimizing  the  overall  absolute  acoustic 
source  strength  of  the  control  sources.  Mathematically,  this  amounts  to  the  minimization  of  multi- variable 
complex-valued  functions  in  the  sense  of  Li  with  conical  constraints,  which  are  only  “marginally”  convex. 
The  corresponding  numerical  optimization  problem  appears  very  challenging  even  for  the  most  sophisticated 
state-of-the-art  methodologies,  and  even  when  the  dimension  of  the  grid  is  small,  and  the  waves  are  long. 

Our  central  result  is  that  the  global  Li-optimal  solution  can,  in  fact,  be  obtained  without  solving  the 
numerical  optimization  problem.  This  solution  is  given  by  a  special  layer  of  monopole  sources  on  the 
perimeter  of  the  protected  region.  We  provide  a  rigorous  proof  of  the  global  Li  minimality  for  both  continuous 
and  discrete  optimization  problems  in  the  one-dimensional  case.  We  also  provide  numerical  evidence  that 
corroborates  our  result  in  the  two-dimensional  case,  when  the  protected  domain  is  a  cylinder.  Even  though 
we  cannot  fully  justify  it,  we  believe  that  the  same  result  holds  in  the  general  case,  i.e.,  for  multi-dimensional 
settings  and  domains  of  arbitrary  shape.  We  formulate  it  as  a  conjecture  at  the  end  of  the  paper. 

Key  words,  noise  cancellation,  control  sources,  minimization  of  amplitude,  surface  monopoles 

Subject  classification.  Applied  and  Numerical  Mathematics 

1.  Introduction.  The  area  of  active  control  of  sound  has  a  rich  history  of  development,  both  as  a 
chapter  of  theoretical  acoustics,  and  in  the  perspective  of  many  different  applications.  Any  attempt  to 
adequately  overview  this  extensive  area  in  the  framework  of  a  focused  research  publication  would  obviously 
be  deficient.  As  such,  we  simply  refer  the  reader  to  monographs  [3,5,9]  that,  among  other  things,  contain  a 
detailed  survey  of  the  literature. 

The  formulation  of  the  problem  that  we  use  in  the  current  paper  has  been  introduced  and  studied  in 
our  previous  work  [7];  here,  we  analyze  this  formulation  from  the  standpoint  of  optimization.  Let  be  a 
given  domain,  C  ,  where  the  space  dimension  n  =  2  or  n  =  3  (these  two  cases  are  most  interesting  for 
applications).  The  domain  Q  can  be  either  bounded  or  unbounded;  for  the  reason  of  simplicity  only  we  will 
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further  assume  that  rt  is  bounded.  Let  F  be  the  boundary  of  F  =  d^.  Both  on  rt  and  on  its  (unbounded) 
complement  Qi  =  we  consider  the  time-harmonic  acoustic  field  u  =  u{x),  x  e  V' ,  governed  by  the 

non-homogeneous  Helmholtz  equation: 


Lu  =  Au  -h  k^u  =  /.  (1.1) 

Equation  (1.1)  is  subject  to  the  Sommerfeld  radiation  boundary  conditions  at  infinity,  which  for  n  =  2  are 
formulated  as 

u{x)  =  O  , 

and  for  n  =  3  as 

u{x)  =  O  (|a:|“^) , 

The  Sommerfeld  boundary  conditions  specify  the  direction  of  wave  propagation,  and  distinguish  between 
the  incoming  and  outgoing  waves  at  infinity  by  prescribing  the  outgoing  direction  only;  they  guarantee  the 
unique  solvability  of  the  Helmholtz  equation  (1.1)  for  any  compactly  supported  right-hand  side  /  =  f{x). 
We  denote  supp/  =  {x\f{x)  ^  0}. 

The  source  terms  /  =  f{x)  in  equation  (1.1)  can  be  located  on  both  Q  and  its  complement 
to  emphasize  the  distinction,  we  denote 


as  \x\  — >  00, 


(1.2a) 


du{x) 

d\x\ 


+  iku(x)  =  o  (|a:| 


as  \x\  — >  00. 


(1.2b) 


/  =  /+  +  /-,  (1.3) 

where  the  sources  /+  are  interior,  supp  /+  C  fl,  and  the  sources  f~  are  exterior,  supp  f~  C  fli,  with  respect 
to  Q.  Accordingly,  the  overall  acoustic  field  u  =  u{x)  can  be  represented  as  the  sum  of  two  components: 

u  =  u~,  (1.4) 

where 

Lu^  =  /+,  (1.5a) 

Lu-  =f-.  (1.5b) 

Note,  both  =  u'^{x)  and  u-  =  u-{x)  are  defined  on  the  entire  W,  the  superscripts  “-h”  and  ”  refer 
to  the  sources  that  drive  each  of  the  field  components  rather  than  to  the  domains  of  these  components.  The 
setup  described  above  is  schematically  shown  in  Figure  1.1. 

Hereafter,  we  will  call  the  component  of  (1.4),  (1.5a)  sound,  or  “friendly”  part  of  the  total  acoustic 
field;  the  component  u-  of  (1.4),  (1.5b)  will  accordingly  be  called  noise,  or  “adverse”  part  of  the  total  acoustic 
field.  In  the  formulation  that  we  are  presenting,  ft  will  be  a  (predetermined)  region  of  space  to  be  shielded. 
This  means  that  we  would  like  to  eliminate  the  noise  inside  Q  while  leaving  the  sound  component  there 
unaltered.  In  the  mathematical  framework  that  we  have  adopted,  the  component  u-  of  the  total  acoustic 
field,  i.e.,  the  response  to  the  adverse  sources  /“  (see  (1.3),  (1.4),  (1.5)),  will  have  to  be  cancelled  on  fl, 
whereas  the  component  W^,  i.e.,  the  response  to  the  friendly  sources  /+,  will  have  to  be  left  unaffected  on 
fl.  A  physically  more  involved  but  conceptually  easy  to  understand  example  that  can  be  given  to  illustrate 
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Fig.  1.1.  Geometric  setup. 


the  foregoing  idea  of  shielding,  is  that  inside  the  passenger  compartment  of  an  aircraft  we  would  like  to 
eliminate  the  noise  coming  from  the  propulsion  system  located  outside  the  fuselage  while  not  interfering 
with  the  ability  of  the  passengers  to  listen  to  the  inflight  entertainment  programs  or  to  converse. 

The  concept  of  active  noise  control  that  we  will  be  discussing  implies  that  the  component  u~  is  to  be 
suppressed  on  fl  by  introducing  additional  sources  of  sound  g  =  g{x)  exterior  with  respect  to  ft,  supp^  C 
so  that  the  total  acoustic  field  u  =  u{x)  be  now  governed  by  the  equation  (cf.  formulae  (1.1),  (1.3)) 

Lu  =  f+  +  f-+g,  (1.6) 


and  coincide  with  only  the  friendly  component  on  the  domain  ft: 


u 


=  u ' 


xe^ 


(1.7) 


The  new  sources  g  =  g{x)  of  (1.6),  see  Figure  1.1,  will  hereafter  be  referred  to  as  the  control  sources 
or  simply  controls.  An  obvious  solution  for  these  control  sources  is  ^  =  —f~.  This  solution,  however,  is 
excessively  expensive.  One  one  hand,  the  excessiveness  comes  from  the  information-type  considerations,  as 
the  solution  g  =  —f~  requires  explicit  and  detailed  knowledge  of  the  structure  and  location  of  the  sources 
f~.  As  shown  in  [7],  this  knowledge  is,  in  fact,  superfluous.  On  the  other  hand,  the  implementation  of  the 
solution  g  —  —f~  may  encounter  most  serious  technical  difficulties.  In  the  previous  example,  it  is  obviously 
not  feasible  to  directly  counter  the  actual  noise  sources,  which  are  aircraft  propellers  or  turbofan  jet  engines 
located  on,  or  underneath,  the  wings.  Therefore,  other  solutions  of  the  foregoing  noise  control  problem. 
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besides  the  most  obvious  one,  may  be  preferable  from  both  the  theoretical  and  practical  standpoint.  The 
general  solution  for  the  control  sources  g  was  obtained  in  our  previous  work  [7],  and  we  describe  it  in  the 
next  Section  2. 

Before  proceeding,  let  us  only  note  that  in  the  current  paper  we  focus  on  the  case  of  the  standard 
constant-coefficient  Helmholtz  equation  (1.1),  which  governs  the  acoustic  field  and  is  valid  throughout  the 
entire  space  W.  This  allows  us  to  make  the  forthcoming  analysis  most  straightforward.  As  a  matter  of 
fact,  other,  more  complex,  cases  that  involve  variable  coefficients  and,  maybe,  nonlinearities  in  the  governing 
equations  over  some  regions,  as  well  as  different  types  of  far-field  behavior,  discontinuities  in  the  material 
properties,  etc.,  can  be  considered  as  well.  Approaches  to  obtaining  solutions  for  active  controls  in  these  cases 
are  outlined  in  our  previous  paper  [7]  for  the  continuous  formulation  of  the  problem,  and  in  the  monograph 
by  Ryaben’kii  [13,  Part  VIII]  for  the  discrete  formulation  of  the  problem. 

The  material  in  the  rest  of  the  paper  is  organized  as  follows.  In  Section  2,  we  introduce  the  control 
sources  for  the  continuous  formulation  of  the  problem.  In  Section  3,  we  obtain  the  control  sources  in  the 
discrete  formulation  of  the  problem,  i.e.,  on  the  grid.  In  Section  4,  we  discuss  minimization  of  the  overall 
acoustic  source  strength  of  active  controls  that  we  have  constructed,  which  mathematically  amounts  to  the 
optimization  in  the  sense  of  Li.  We  present  a  convincing  two-dimensional  numerical  evidence,  as  well  as  a 
rigorous  one-dimensional  proof,  of  global  Li  optimality  of  a  particular  layer  of  monopole  sources  concentrated 
on  the  perimeter  of  the  protected  region.  We  believe  that  the  combination  of  computations  in  two  space 
dimensions  and  general  proof  in  one  space  dimension  cannot  be  a  mere  coincidence.  As  such,  even  though 
we  cannot  fully  justify  it,  we  formulate  the  corresponding  general  result  on  global  Li  optimality  of  surface 
monopoles  as  a  conjecture  in  the  concluding  Section  5.  It  basically  implies  that  the  aforementioned  Li 
optimization  problem  can  be  solved  without  using  any  numerical  optimization  techniques. 

2.  Continuous  Control  Sources. 


2.1.  General  Solution.  As  demonstrated  in  [7],  the  general  solution  for  the  control  sources  g  =  g{x) 
is  given  by  the  following  formula 


g{x)  =  -Lw 


(2.1) 


where  w  =  w{x),  x  G  f7i,isa  special  auxiliary  function-parameter  that  parameterizes  the  family  of  controls 
(2.1).  The  requirements  that  the  function  w{x)  must  meet  are,  in  fact,  relatively  “loose.”  At  infinity,  it  has 
to  satisfy  the  Sommerfeld  boundary  conditions  (1.2a)  or  (1.2b).  At  the  interface  T,  the  function  w  and  its 
normal  derivative  have  to  coincide  with  the  corresponding  quantities  that  pertain  to  the  total  acoustic  field 
u  given  by  formula  (1.4)  d 


dw 

du 

=  u 
r 

r’  dn 

P  dn 

(2.2) 


Other  than  that,  the  function  w{x)  used  in  (2.1)  is  arbitrary,  and  consequently  formula  (2.1)  defines  a  large 
family  of  control  sources,  which,  as  will  be  seen  in  Section  4,  provides  ample  room  for  optimization.  To  make 
the  discussion  in  the  current  paper  self-contained,  we  briefiy  outline  below  the  justification  for  formula  (2.1) 
as  general  solution  for  controls,  while  referring  the  reader  to  [7]  for  further  detail. 


^In  practice,  the  quantities  u  and  ^  on  F  can  be  measured  and  supplied 
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Introducing  the  fundamental  solution  G  =  G{x)  of  the  Helmholtz  operator  L  of  (1.1)  for  n  =  2: 

G{x)  =  -^.4^\k\x\),  (2.3a) 

where  =  Jo{z)  —  iYo{z)  is  the  Hankel  function  of  the  second  kind  ,  and  for  n  =  3: 

p-ik\x\ 

we  can  obviously  represent  the  sound  portion  =  u'^{x)  of  the  overall  acoustic  field  that  satisfies  equation 
(1.5a)  everywhere  on  as  follows: 

a+(a:)  =  [  G{x  -  y)f+{y)dy  =  [  G{x  -  y)Lu^{y)dy,  a:  e  M".  (2.4) 

Consequently,  applying  the  classical  Green’s  formula  (see,  e.g.,  [16]  or  [22])  to  the  function  u'^  =  u'^{x)  on 
we  have: 


dn 


:  0,  xen, 


(2.5) 


where  integrals  in  (2.5)  are,  again,  convolutions.  Similarly,  applying  the  same  Green’s  formula  on  to 
u~  =  u~{x)  and  using  (2.5),  we  obtain: 


u  {x)-J^(u  Iri'^dn  dn^)  ®  ^ 

Therefore,  from  (2.6)  we  can  conclude  that  the  desired  annihilating  acoustic  signal  v  =  v{x)  that  cancels 
out  the  unwanted  noise  on  can  be  obtained  as 

so  that  foi  X  e  ft  we  indeed  have 


\u{x)-\-v{x)  =u'^{x) 

L  J  xgQ 


To  actually  implement  the  annihilating  signal  v  of  (2.7),  we  introduce  the  auxiliary  function  w  =  w{x) 
on  W  that  satisfies  the  aforementioned  conditions  at  the  interface  T  and  at  infinity,  and  apply  the  Green’s 
formula  to  which  yields: 


(x)  —  /  GLwdy 
Jq 


=  /r(“’ 


dG  dw 
dn  dn 


dsy,  X  eft. 


(2.8) 


As  w{x)  satisfies  the  Sommerfeld  conditions  (1.2a)  or  (1.2b),  we  obviously  have  w{x)  =  GLwdy ^  x 

which,  along  with  formulae  (2.2)  and  (2.7),  allows  us  to  transform  equality  (2.8)  to 


v{x)  =  —  f  GLwdy.  (2.9) 

J 

Equation  (2.9)  implies  that  for  any  w  =  w{x)  chosen  as  described  above,  formula  (2.1)  describes  an  appro¬ 
priate  control  function. 
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Conversely,  assume  that  g  =  9{x)^  supp^  C  is  a  control  field  such  that  the  solution  u  =  u{x) 
of  equation  (1.6)  subject  to  the  Sommerfeld  conditions  (1.2a)  or  (1.2b)  satisfies  equality  (1.7).  Then,  by 
choosing  =  ^^;(^c),  a;  G  as  the  solution  to  the  nonhomogeneous  equation 

-Lw  =  g- 

subject  to  the  corresponding  Sommerfeld  condition,  (1.2a)  or  (1.2b),  one  can  represent  g  in  the  form  (2.1), 
see  [7].  Altogether  we  obtain  that  formula  (2.1)  describes  the  general  solution  for  controls.  In  other  words, 
for  any  ^^;(^c),  formula  (2.1)  provides  an  appropriate  control  field  g{x)^  and  any  appropriate  control  field  g{x) 
can  be  represented  in  the  form  (2.1)  with  some  particular  choice  of  w{x). 

Let  us  emphasize  several  important  properties  of  controls  (2.1).  First  of  all,  from  the  foregoing  derivation 
we  can  see  that  to  obtain  these  controls  one  needs  no  knowledge  of  the  actual  exterior  sources  of  noise  /“. 
In  other  words,  neither  their  location,  nor  structure,  nor  strength  are  required.  All  one  needs  to  know  is 
u  and  on  the  perimeter  F  of  the  protected  region  fl.  As  has  been  mentioned,  in  a  practical  setting 

and  can  be  interpreted  as  measurable  quantities  that  are  supplied  to  the  control  system  as  the 
input  data.  Moreover,  these  measurable  quantities  refer  to  the  overall  acoustic  field  u,  rather  than  only  to 
its  unwanted  component  u~.  In  other  words,  the  methodology  can  automatically  distinguish  between  the 
signals  coming  from  the  exterior  and  interior  sources,  and  can  tune  the  controls  so  that  they  cancel  only 
the  unwanted  exterior  signal.  This  capability  is  extremely  important  as  in  many  applications  the  overall 
acoustic  field  always  contains  a  component  that  needs  to  be  suppressed  along  with  the  part  that  needs  to 
be  left  intact.  Alternatively,  one  can  say  that  the  control  sources  (2.1)  are  insensitive  to  the  interior  sound 
u'^{x).  Indeed,  given  a  function  w{x)  that  satisfies  interface  conditions  (2.2)  and  the  radiation  boundary 
conditions  at  infinity,  we  can  take  instead  w{x)  =  w{x)  —  u'^{x)]  this  new  function  will  satisfy  the  interface 
conditions  (2.2)  with  u  replaced  by  u~  and  the  same  Sommerfeld  conditions  at  infinity.  Most  important, 
the  control  sources  generated  by  w{x)  will  be  the  exact  same  control  sources  as  those  generated  by  w{x) 
because  Lw{x)  =  L[w{x)  —  u'^{x)]  =  Lw{x)  for  x  e  fti. 

Let  us  also  note  that  in  a  more  general  framework,  formula  (2.9)  (with  (2.2)  taken  into  account)  can 
be  interpreted  as  a  particular  case  of  the  generalized  potential  of  Calderon’s  type  with  the  vector  density 
—  (ii,  1^)  .  This  more  general  framework  allows  us  to  analyze  more  complex  formulations  of  the  active 

noise  control  problem,  see  [7],  such  as  those  that  involve  variations  in  material  properties  and  alternative 
types  of  far-field  behavior  of  the  solution.  We  refer  the  reader  to  the  original  work  by  Calderon  [2]  and 
Seeley  [14],  as  well  as  to  the  monograph  by  Ryaben’kii  [13],  for  the  general  concepts  related  to  Calderon’s 
potentials  and  associated  pseudodifferential  boundary  projection  operators.  As  concerns  the  aforementioned 
more  advanced  formulations  of  the  noise  control  problem,  the  key  result  is  basically  the  same  as  above. 
The  general  solution  for  control  sources  is  still  given  by  formula  (2.1),  where  the  auxiliary  function  w{x) 
should  still  satisfy  the  interface  conditions  (2.2)  and  the  problem-specific  far-field  boundary  conditions  (in 
case  they  differ  from  the  previously  mentioned  Sommerfeld  conditions).  The  operator  L  in  formula  (2.1)  will, 
however,  no  longer  be  the  constant-coefficient  Helmholtz  operator  of  (1.1),  it  will  rather  be  the  problem- 
specific  variable-coefficients  operator  that  accounts  for  the  particular  variations  in  material  properties,  etc. 
As,  however,  formula  (2.1)  for  controls  does  not  change  (see  [7]),  we  immediately  conclude  that  we,  in  fact, 
do  not  need  to  know  the  operator  L  on  fl.  In  other  words,  for  obtaining  the  control  sources  we  do  not  need 
to  know  the  material  properties  of  the  sound-conducting  medium  inside  the  protected  region.  This  result. 
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see  [7],  which  first  seems  counterintuitive,  has,  in  fact,  a  straightforward  physical  explanation.  One  only 
needs  to  realize  that  the  noise  we  want  to  suppress,  and  the  output  of  controls  that  is  supposed  to  annihilate 
this  noise,  propagate  across  one  and  the  same  medium,  and  we  do  not  need  to  know  what  this  medium  is 
(under  only  some  relatively  non-rest rictive  limitations,  see  [7]).  It  is  also  interesting  to  mention  that  the 
aforementioned  Calderon’s  boundary  projection  operators  essentially  render  the  decomposition  of  the  wave 
field  u{x)  on  the  boundary  F  into  its  incoming  and  outgoing  component  with  respect  to  the  domain  ft, 
see  [7].  Subsequently,  the  controls  (2.1)  can  be  interpreted  as  either  sources  cancelling  the  incoming  wave 
field  for  the  domain  to  be  shielded,  i.e.,  ft,  or  alternatively,  as  sources  cancelling  the  outgoing  wave  field  for 
the  domain  complementary  to  the  one  to  be  shielded,  i.e.,  fti  =  The  latter  interpretation  is  often 

more  versatile,  see  [7]. 

Another  important  thing  to  notice  is  that  the  control  sources  g{x)  of  (2.1)  are  defined,  generally  speaking, 
on  the  entire  complementary  domain  fti  =  For  the  analysis  of  specific  problems,  especially  when 

the  protected  region  ft  is  bounded  and  the  complementary  region  fti  is  unbounded,  as  in  Figure  1.1,  it 
may  be  convenient  to  consider  compactly  supported  control  sources,  i.e.,  the  control  sources  concentrated 
in  the  vicinity  of  the  interface  F.  To  obtain  such  controls,  one  needs  to  narrow  down  the  class  of  functions 
w{x)  used  in  formula  (2.1).  Namely,  instead  of  considering  arbitrary  w{x)  subject  only  to  constraints  (2.2) 
and  Sommerfeld  boundary  conditions  at  infinity,  one  needs  to  consider  w{x)  that  become  a  solution  to  the 
homogeneous  Helmholtz  equation  everywhere  outside  some  larger  domain  that  fully  contains  the  protected 
region  ft.  In  so  doing,  the  area  outside  ft  that  supports  the  controls  g{x)  of  (2.1)  will  stretch  from  F  to  the 
outer  boundary  of  the  aforementioned  larger  domain,  and  may  basically  look  as  a  curvilinear  strip  adjacent 
to  F  from  the  exterior  side.  This  strip  may,  in  principle,  be  made  as  narrow  as  desired  and  may  eventually 
shrink  completely  thus  reducing  only  to  the  interface  F  itself.  As  will  be  seen,  the  sources  g{x)  supported 
only  on  the  perimeter  F,  represent  an  important  class  of  active  controls.  Such  distributions  g{x)  include 
monopole  and  dipole  layers  as  special  cases,  which  are  idealizations  of  pulsating  or  vibrating  membranes. 

2.2.  Artificial  Boundary  Conditions  and  Compactly  Supported  Controls.  To  actually  obtain 
compactly  supported  controls  g{x)  in  a  particular  setting,  it  is  often  convenient  to  use  the  methodology 
known  as  artificial  boundary  conditions  (see  the  review  paper  [17])  for  the  selection  of  the  appropriate 
function  w{x)  in  formula  (2.1).  Assume  that  w{x)  satisfies  the  homogeneous  Helmholtz  equation  Lw  =  0 
outside  some  outer  artificial  boundary,  which  is  a  closed  surface  (curve)  with  the  interior  that  fully  contains 
ft.  In  Figure  1.1,  we  schematically  represent  this  outer  boundary  as  a  sphere  (circle)  of  radius  R.  It  turns  out 
that  one  can  equivalently  replace  the  homogeneous  equation  Lw  =  0  along  with  the  Sommerfeld  boundary 
conditions  at  infinity  by  the  special  artificial  boundary  conditions  (ABCs)  at  the  outer  boundary.  For  the 
outer  boundary  of  a  general  shape,  this  can  be  done  most  efficiently  using  the  same  apparatus  of  Calderon’s 
pseudodifferential  boundary  projection  operators,  see  [13,17],  that  has  been  mentioned  before.  For  the 
particular  case  of  a  regular  spherical  or  circular  outer  boundary,  see  Figure  1.1,  which  is  convenient  because 
of  the  simplicity  of  the  forthcoming  analysis,  the  ABCs  can  be  constructed  as  follows. 

We  will  use  spherical  coordinates  (y9,0, 0)  in  n  =  3,  and  assume  that  Lw  =  0  foi  \x\  =  p  >  R, 
see  Figure  1.1.  In  addition,  we  will  assume  that  w{x)  satisfies  the  Sommerfeld  boundary  condition  (1.2b). 
Expanding  w{x)  with  respect  to  spherical  functions  Yi^,  /  =  0, 1, 2, . . . ,  m  =  0, ±1, ...  ,±l,  and  separating 
the  variables  in  the  differential  operator  L,  we  arrive  at  the  following  collection  of  the  second-order  ordinary 
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differential  equations 


2  dWifYi  2 

~d^  p  [ 

p  ^  Rf  /  =  0, 1, 2, . . .  , 


/(/  +  !)' 
p^  . 


'dJlm  —  0? 


m  =  0,  ±1, . . .  ,  zb/, 


(2.10) 


for  the  unknown  radial  modes  wim-  These  modes  are  also  supposed  to  satisfy  boundary  conditions  at  infinity 

Wim{p)  =  o  {p~^) ,  +  ikwimip)  =  O  {p~^) ,  asp — ^  (X),  (2.11) 

which  immediately  follow  from  the  Sommerfeld  condition  (1.2b).  For  any  given  pair  (/,m)  the  general 
solution  of  equation  (2.10)  is  given  by 

Wlm  =  (2.12) 

where  ci  and  C2  are  arbitrary  constants,  and  and  are  the  Hankel  functions  of  the  first 

and  second  kind,  respectively.  Taking  into  account  the  asymptotic  expressions  for  the  Hankel  functions  for 
a  fixed  order  u  and  large  values  of  p  (see,  e.g.,  [22]): 


K''  -  5*'  -  91 + °  ' 

=  \/5“p  H  (p  “  f"  ■  91  +  °  ■ 

we  conclude  that  only  p~^^‘^Hi^-^i^{kp)  satisfies  boundary  conditions  (2.11)  and  consequently,  the  constant 
Cl  (see  (2.12))  in  any  particular  solution  that  satisfies  the  Sommerfeld  conditions  at  infinity  has  to  be  equal 
to  zero.  As  the  two  functions  p~^^‘^H^^^i^{kp)  and  p~^^‘^H^^^i^{kp)  form  a  fundamental  system  of  solutions 
for  the  linear  homogeneous  second-order  ODE  (2.10),  requiring  that  only  one  of  them,  p~^^^H^^^i^{kp)^  be 
present  in  the  actual  solution  wim  is  equivalent  to  requiring  that  wimip)  be  parallel  to  for 

p>  Rin  the  sense  of  the  corresponding  Wronskian  vanishing  at  p  =  R: 


det 


'dJlm 


(2.13) 


Obviously,  equality  (2.13)  enforced  at  y9  =  i?  implies  that  it  will  hold  for  all  y9  >  i?  as  well.  Equality  (2.13) 
is  a  linear  homogeneous  relation  between  a  given  Fourier  component  wim  of  the  solution  w{x)  and  the 
corresponding  Fourier  component  -^wim  of  its  normal  derivative  -^w{x)  on  the  spherical  surface  p  =  R. 
The  entire  family  of  such  relations  for  all  /  =  0, 1, 2, . . .  and  m  =  0,  ±1, . . .  ,  ±/  set  at  y9  =  i?  is  equivalent 
to  saying  that  the  function  w  =  w{x)  originally  defined  inside  the  sphere,  i.e.,  for  p  <  R,  can  be  smoothly 
extended  to  the  region  p  >  R  so  that  the  extension  will  solve  the  homogeneous  equation  Lw  =  0  foi  p>  R 
and  have  a  proper  far-held  behavior,  i.e,  satisfy  the  Sommerfeld  condition  (1.2b).  Hereafter,  we  will  refer 
to  relations  (2.13)  for  /  =  0, 1, 2, . . .  and  m  =  0,  zbl, . . .  ,  zb/  as  to  the  artificial  boundary  conditions  for  the 
three-dimensional  Helmholtz  equation  on  the  spherical  surface  p  =  R. 

The  ABCs  for  the  two-dimensional  case  on  the  circular  external  artihcial  boundary  p  =  R  can  be 
obtained  similarly.  We  introduce  polar  coordinates  (p^O)  in  n  =  2,  use  standard  Fourier  expansion  in 
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the  circumferential  direction  with  respect  to  the  complex  exponents  e  =  0,  ±1,  ±2, . . . ,  and  arrive  at 
the  following  collection  of  the  second-order  ordinary  differential  equations: 


1  dwi 

- 1  - 1  + 

dp^  p  dp 


wi  =  0, 


p>  R,  1  =  0,  ±1,  ±2, . . .  , 


(2.14) 


for  the  unknown  radial  modes  wi .  These  modes  are  also  supposed  to  satisfy  boundary  conditions  at  infinity 

Wi{p)  =  O  +  ikwi{p)  =  o  ,  as  /9  — ^00,  (2.15) 

which  immediately  follow  from  the  Sommerfeld  condition  (1.2a).  For  every  given  I,  equation  (2.14)  is  the 
Bessel  equation  and  has  general  solution 


Wi  =  ciH^^\kp)  +  C2Hl^\kp), 


(2.16) 


where  ci  and  C2  are,  again,  arbitrary  constants.  The  asymptotics  of  the  Hankel  functions  for  large  p^s 
indicates  that  to  satisfy  (2.15)  one  must  have  ci  =  0  in  any  particular  solution  that  satisfies  the  radiation 
conditions  at  infinity.  This  requirement  leads  to  the  following  ABCs  for  the  two-dimensional  Helmholtz 
equation: 


wi  (kp) 


(2.17) 


where  relations  (2.17)  should  be  considered  for  all  I  =  0,  ±l,zb2, -  We  refer  the  reader  to  the  review 

article  [17]  for  further  detail  on  the  construction  of  ABCs  for  different  equations  in  different  settings.  Let 
us  also  reiterate  that  once  the  ABCs  (2.13)  or  (2.17)  are  satisfied  for  all  radial  modes,  then  we  can  consider 
Lw  =  0  for  y9  >  i?,  and  as  such,  the  resulting  control  sources  g{x)  given  by  (2.1)  will  be  compactly  supported 
between  the  interface  F  and  the  external  artificial  boundary  p  =  R. 


2.3.  Types  of  Control  Sources.  Let  us  now  analyze  the  continuous  control  sources  from  the  stand¬ 
point  of  their  geometric  location  and  type  of  acoustic  excitation  that  they  provide.  To  put  this  analysis  into 
a  mathematical  perspective,  it  will  be  convenient  to  use  the  apparatus  of  distributions,  see,  e.g.,  [22]. 

In  our  original  derivation  of  formula  (2.1),  we  have  implicitly  assumed  that  the  function  w{x)  was 
sufficiently  smooth  so  that  the  operator  L  could  be  applied  in  the  classical  sense  everywhere  on  fti  = 

In  this  case,  the  function  g{x)  is  locally  integrable,  and  can  be  interpreted  as  regular  distribution.  As 
any  other  distribution,  it  can  be  represented  as  convolution  of  its  own  self  with  the  5-function:  g  =  S  ^  g, 
which  means  that  the  control  field  g(x)  can  be  formally  viewed  as  the  “sum”  of  elementary  point  monopoles 
g{y)S{x  —  y)  with  regular  density  g:  g{x)  =  5^g  =  J  g{y)S{x  —  y)dy.  Next,  we  remind  that  by  definition  of 
the  fundamental  solution:  LG  =  5,  the  response  to  every  such  elementary  monopole  g{y)^{x  —  y)  will  be  given 
9{y)G{x  —  y)]  and  consequently,  the  overall  control  output  G^g  =  J g{y)G{x  —  y)dy  shall  be  interpreted 
as  superposition  of  the  foregoing  elementary  responses,  i.e.,  solutions  generated  by  the  aforementioned  point 
monopoles.  Altogether  we  see  that  the  original  formula  (2.1)  provides  the  general  solution  for  controls  in  a 
particular  class  of  functions,  namely,  that  of  regular  distributions,  which  in  terms  of  physics  is  characterized 
as  volumetric  control  sources  of  the  monopole  type  on  the  complementary  domain  Compactly 

supported  control  sources  discussed  in  Section  2.2  obviously  fall  into  this  category. 
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In  many  cases  it  may  also  be  desirable  to  consider  surface  controls,  i.e.,  the  control  sources  that  are 
concentrated  only  on  the  interface  F.  Let  us  first  assume  that  there  are  no  interior  sources,  which  means  that 
u'^(x)  =  0,  and  the  acoustic  field  we  want  to  control  consists  only  of  its  adverse  component,  u{x)  =  u~{x). 
After  the  control,  it  has  to  be  equal  to  zero  on  the  domain  ft.  As  shown  in  [18],  the  general  solution  for 
surface  controls  is  given  by 

^^(0- J^(k-«H(r)),  (2.18) 

where  w  =  as  before,  denotes  the  auxiliary  function-parameter  that  in  this  case  has  to  satisfy  the 

homogeneous  Helmholtz  equation  on  the  complementary  domain:  Lw  =  0  for  cc  G  ili,  and  the  Sommerfeld 
boundary  condition  (1.2a)  or  (1.2b)  at  infinity.  Expressions  in  rectangular  brackets  in  formula  (2.18)  denote 
discontinuities  of  the  corresponding  quantities  across  the  interface  F.  In  so  doing,  the  first  term  on  the 
right-hand  side  of  (2.18)  represents  the  density  of  a  single-layer  potential,  which  is  a  layer  of  monopoles  on 
the  interface  F,  and  the  second  term  on  the  right-hand  side  of  (2.18)  represents  the  density  of  a  double-layer 
potential,  which  is  a  layer  of  dipoles  on  the  interface  F. 

A  detailed  justification  of  formula  (2.18)  as  general  solution  for  surface  controls  can  be  found  in  [18]. 
Here  we  only  mention  that  it  basically  amounts  to  proving  that  a  given  solution  u{x)  of  the  homogeneous 
equation  Lu  =  0  on  Q  can  be  represented  as  a  combination  of  a  single-layer  potential  and  a  double-layer 
potential  if  and  only  if  the  densities  of  the  aforementioned  potentials  are  defined  as  —  |^]j^5(F)  and 
^  {[w  —  ii]r5(F)),  respectively  (cf.  formula  (2.18)).  The  direct  implication  is  easy  to  establish  by  applying 
the  operator  L  to  the  discontinuous  function 


(surf)  _  _ 


dw  du 
dn  dn 


v{x) 


u{x)  for  X 
w{x)  for  X  e  fli 


in  the  sense  of  distributions,  see  [22],  which  yields 


Lv  = 


dw 

dn 


(2.19) 


(2.20) 


and  subsequently  reconstructing  M(a:)|  as  convolution  with  the  fundamental  solution.  The 

inverse  implication  requires  explicitly  obtaining  ^^;(^c)  for  a  given  u{x)  and  given  surface  densities,  which  is 
done  in  [18],  again,  in  the  form  of  a  special  surface  integral.  Then,  the  control  sources  (2.18)  are  obtained 
by  simply  taking  (2.20)  with  the  opposite  sign,  which  guarantees  the  cancellation  of  u{x)  on  Q  by  —v{x). 

In  the  family  of  surface  controls  (2.18)  we  identify  two  important  particular  cases.  First,  the  cancellation 
of  u{x),  X  ^  can  be  achieved  by  using  only  the  surface  monopoles,  i.e.,  by  employing  only  a  single-layer 
potential  as  the  annihilating  signal.  To  do  that,  we  need  to  find  ^^;(^c),  a;  G  ili,  such  that  the  overall  function 
v{x)  of  (2.19)  would  have  the  discontinuity  on  F  only  in  its  normal  derivative,  and  not  in  the  function  itself. 
This  w{x)  will  then  be  a  solution  of  the  following  external  Dirichlet  problem: 


Lw  =  0,  X 


(2.21) 


w\ 


=  u 


r’ 


subject  to  the  appropriate  Sommerfeld  boundary  condition  (1.2a)  or  (1.2b).  Problem  (2.21)  is  always 
uniquely  solvable  on  fti  =  Second,  one  can  employ  only  the  double-layer  potential  to  cancel  out 
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u{x),  X  ^  i.e.,  use  only  surface  dipoles  as  the  control  sources.  In  this  case,  the  function  x  G 

has  to  be  such  that  v{x)  given  by  (2.19)  would  have  discontinuity  on  F  only  in  the  function  itself,  and  not 
in  its  normal  derivative.  This  w{x)  should  then  solve  the  following  external  Neumann  problem: 

Lw  =  0,  a;  G 

dw  _  du  (2.22) 

dn  r  dn  r’ 

again,  subject  to  the  appropriate  Sommerfeld  condition  at  infinity,  (1.2a)  or  (1.2b);  the  latter  guarantees  the 
solvability  of  (2.22). 

In  a  more  general  case,  when  interior  sources  are  present,  the  results,  in  fact,  do  not  change.  Assume, 
as  before,  that  the  overall  acoustic  field  is  the  sum  of  its  friendly  and  adverse  components  (see  (1.4)): 
u{x)  =  u^{x)  +  u~{x)^  and  we  want  to  cancel  out  u~{x)  on  fl.  Then,  formula  (2.18),  where  u  shall  now 
be  interpreted  as  in  (1.4)  and  is  a  function-parameter,  will  still  provide  the  general  solution  for  surface 
controls.  Indeed,  since  Lu'^  =  0  on  and  moreover,  ii+(x)  and  ^^(ic)  are  continuous  across  the  interface 
r,  then  (2.18)  simply  reduces  to 

=  -  [I!  -  fir] , 

where  the  new  function-parameter  w  is  given  by  w{x)  =  w{x)  —  u'^{x)  and  as  such,  satisfies  the  afore¬ 
mentioned  general  requirements  of  ^^;’s.  This  means  that  the  control  sources  cc  G  F,  defined  by 

(2.18),  or  equivalently  (2.23),  appear  insensitive  to  the  friendly  component  u^{x)  of  the  acoustic  field,  and 
will  annihilate  precisely  u~{x)  on  the  domain  Vt.  Furthermore,  the  entire  family  of  surface  controls  (2.23)  is 
obviously  the  exact  same  family  as  we  would  have  obtained  if  there  were  no  interior  sources  and  the  overall 
acoustic  field  consisted  of  only  u~{x).  It  is  also  clear,  that  the  same  reasoning  will  apply  to  the  particular 
cases  of  purely  monopole  and  purely  dipole  controls.  Namely,  if  in  problems  (2.21)  and  (2.22)  we  interpreted 
the  boundary  data  u  and  (respectively)  in  the  sense  of  (1.4),  then  the  solution  w  of  either  problem  would 
obviously  be^^;(^c)  =  w{x)  u'^{x),  where  w{x)  is  the  solution  that  corresponds  to  u'^{x)  =0.  This  means 
that  both  the  monopole  and  the  dipole  layer  constructed  using  the  respective  solution  w{x)  would  be  the 
exact  same  monopole  or  dipole  layer  that  suppresses  u~{x)  on  Q  in  the  case  of  no  interior  sources  and  no 
interior  sound. 

Altogether  we  conclude  that  as  indicated  by  formula  (2.18),  surface  control  sources  are  combinations 
of  monopole  and  dipole  layers,  with  the  two  “extreme”  cases  corresponding  to  either  only  monopoles,  see 
(2.21),  or  only  dipoles,  see  (2.22).  From  the  standpoint  of  physics  and  engineering,  the  monopole  and  dipole 
sources  provide  different  types  of  excitation  to  the  surrounding  sound- conducting  medium.  A  point  monopole 
source  can  be  interpreted  as  a  vanishingly  small  pulsating  sphere  that  radiates  acoustic  waves  symmetrically 
in  all  directions,  whereas  a  dipole  source  resembles  more  a  small  oscillating  membrane  that  has  a  particular 
directivity  of  radiation  (see,  e.g.,  [9]  for  a  more  detailed  discussion  on  the  properties  of  different  sources). 

In  the  current  paper,  we  are  going  to  study  the  control  sources  of  one  physical  kind  only;  namely,  the 
monopoles.^  This  class  includes  the  volumetric  controls  (2.1),  i.e.,  monopoles  distributed  in  space,  and  their 

^  We  intentionally  leave  out  the  discussion  of  how  easy  or  how  difficult  it  may  be  to  implement  a  given  type  of  control  sources 
in  a  practical  setting. 
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limiting  case  given  by  the  monopole  layer  on  the  surface  T  (cf.  (2.18)): 


(surf)  _ 

^monopole 


dw  du 
dn  dn 


(2.24) 


where  w  =  w{x)  in  (2.24)  solves  the  Dirichlet  problem  (2.21).  The  reason  for  narrowing  down  the  class  of 
admissible  controls  is  that  having  only  one  and  the  same  physical  type  of  sources  greatly  facilitates  their 
description  and  comparison  using  the  same  quantitative  characteristics,  e.g.,  acoustic  strength  (see  Section  4 
for  more  detail).  This,  in  particular,  will  allow  us  to  formulate  (see  Section  4.1)  and  solve  (see  Sections  4.2 
and  4.3)  the  optimization  problem  for  the  control  sources. 


3.  Discrete  Control  Sources.  Similarly  to  the  continuous  constructions  of  the  previous  Section  2, 
one  can  discretize  the  problem  on  the  grid  and  obtain  the  control  sources  for  the  discrete  formulation.  From 
the  standpoint  of  applications  this  is,  of  course,  preferred,  because  any  practical  design  of  a  noise  control 
system  can  only  contain  a  finite  number  of  elements  or  devices  (acoustic  sensors  and  actuators)  that  will  be 
associated  with  the  grid  nodes  in  the  discrete  case.  Details  regarding  the  discrete  formulation  of  the  noise 
control  problem  can  be  found  in  the  monograph  by  Ryaben’kii  [13,  Part  VIII],  as  well  as  in  the  papers  [20,21]; 
here  we  only  provide  a  brief  account  of  the  corresponding  work.  The  analysis  hereafter  will  not  be  limited  to 
any  specific  type  of  the  grid.  In  particular,  no  adaptation  or  grid  fitting  to  either  the  shape  of  the  protected 
region  (i.e.,  interface  T),  or  that  of  the  external  artificial  boundary,  will  generally  be  required.  In  some 
cases,  though,  it  may  simply  be  convenient  and  inexpensive  to  use  regular  grid  of  an  appropriate  geometry. 
For  example,  as  we  discuss  later,  having  a  polar  or  spherical  grid  may  greatly  simplify  setting  the  ABCs  on 
the  circle  or  sphere  of  radius  i?,  respectively,  in  the  discrete  framework. 


3.1.  Grids  and  Discretization.  Let  us  now  introduce  a  finite-difference  grid  N  that  would  span  both 
fl  and  fli.  In  the  discrete  formulation,  the  grid  never  stretches  all  the  way  to  infinity,  it  is  always  truncated  by 
the  external  artificial  boundary,  which  implies  that  the  discrete  control  sources  that  we  obtain  will  always  be 
compactly  supported.  Later  in  Section  3.3,  we  will  discuss  how  to  set  the  appropriate  ABCs  for  the  discrete 
formulation.  Let  now  be  a  representation  of  the  acoustic  field  on  the  grid,  and  be  a  finite-difference 
approximation  of  the  differential  operator  X  of  (1.1).  To  accurately  define  the  approximation,  we  will  need 
to  introduce  another  grid  M  along  with  the  previously  defined  N.  On  the  grid  M,  we  will  consider  the 
residuals  of  the  operator  ,  and  subsequently  the  right-hand  sides  to  the  corresponding  inhomogeneous 
finite-difference  equation.  We  will  use  the  notations  n  and  m  for  the  individual  nodes  of  the  grids  N  and  M, 
respectively,  and  the  notation  Nm  for  the  stencil  of  the  discrete  operator  X^'^^  centered  at  a  given  m  G  M,  so 
that 


=  E 

nGNm 


(h) 


(3.1) 


where  are  the  coefficients  associated  with  particular  nodes  of  the  stencil.  There  are  no  limitations  to 
the  type  of  discrete  operators  that  one  may  use.  We  only  require  that  the  difference  operator  X^^^  of  (3.1) 
approximate  the  differential  operator  X  of  (1.1)  with  the  accuracy  sufficient  for  a  particular  application. 

Next,  we  introduce  the  following  subsets  of  the  grids  M  and  N,  which  will  allow  us  to  accurately  dis¬ 
tinguish  between  the  interior  and  exterior  domains,  interior  and  exterior  sources,  and  interior  and  exterior 
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solutions  on  the  discrete  level: 

=  M  n  D,  M“  =  M\M+  =  M  n  Di, 

N+  =  U  N-  =  U  (3.2) 

mGMT*'  mGM“ 

7  =  N+  n  N“ ,  7+  =  N“  n  ^7,  7“  =  N+  n 

We  emphasize,  that  the  grid  M  that  pertains  to  the  residuals  of  the  finite-difference  operator  is  parti¬ 
tioned  into  M+  and  M“  directly,  i.e.,  following  the  geometry  of  Q  and  fli.  In  contradistinction  to  that,  the 
grid  N  is  not  partitioned  directly,  we  rather  consider  the  collection  of  all  nodes  of  N  swept  by  the  stencil 
Nrn  when  its  center  belongs  to  M+,  and  call  this  sub-grid  N+,  see  (3.2).  Obviously,  some  of  the  nodes  of 
N+  obtained  by  this  approach  happen  to  be  outside  f),  i.e.,  in  ili,  and  these  nodes  are  called  7“.  The  sets 
N“  and  7+  are  defined  similarly  starting  from  M“ .  The  key  idea  is  that  whereas  the  grids  M+  and  M“ 
do  not  overlap,  the  grids  N+  and  N“  do  overlap,  and  their  overlap  is  denoted  7;  obviously,  7  =  7+07“. 
The  subset  of  grid  nodes  7  is  called  the  grid  boundary^  it  is  a  fringe  of  nodes  that  is  located  near  the  con¬ 
tinuous  boundary  T  and  in  some  sense  straddles  it.  The  specific  structure  of  7  obviously  depends  on  the 
construction  of  the  operator  of  (3.1)  and  the  stencil  Nm  •  For  example,  for  the  conventional  second-order 
central-difference  Laplacians  on  rectangular  grids,  7  will  be  a  two-layer  fringe  of  grid  nodes  located  near  T, 
as  shown  schematically  in  Figure  3.1.  Further  specifics  on  the  construction  of  grid  boundaries  can  be  found 
in  the  monograph  [13]. 

3.2.  Discrete  Noise  Control  Problem  and  Its  General  Solution.  Having  introduced  the  dis¬ 
cretization  (3.1)  and  grid  subsets  (3.2),  we  can  formulate  and  solve  the  noise  control  problem  on  the  grid. 
We  will  reproduce  below  the  key  results  of  [13,  Part  VIII],  see  also  [20,21]  for  detail. 

The  discrete  noise  control  problem  is  formulated  similarly  to  the  continuous  one,  see  Section  1.  Let  /m 
m  G  M+ ,  and  fm^~ ■>  m  G  M“ ,  be  the  interior  and  exterior  discrete  acoustic  sources,  respectively.  Let  ^ 
n  G  N,  and  u"n^~ ,  n  G  N,  be  the  corresponding  solutions,  i.e.,  =  /^^^+  and 

Using  the  same  terminology  as  before,  we  will  call  the  discrete  sound  and  the  discrete  noise. 

The  overall  discrete  acoustic  field  is  the  sum  of  its  sound  and  noise  components,  on 

N,  and  obviously  satisfies  the  equation  =  /^^^+  +  The  goal  is  to  obtain  the  discrete 

control  sources  so  that  the  solution  of  the  equation  =  /^^^+  +  g^^^  be  equal 

to  only  the  sound  component  1^^^^+  on  the  sub-grid  N+ . 

Let  us  now  remind  that  in  the  continuous  case,  the  unique  solvability  of  the  governing  differential 
equation  (inhomogeneous  Helmholtz’  equation)  was  guaranteed  by  the  Sommerfeld  radiation  conditions 
(1.2a)  or  (1.2b)  at  infinity.  In  the  discrete  case,  we  also  need  to  guarantee  the  unique  solvability  of  the 
foregoing  finite-difference  equations,  but  we  obviously  cannot  directly  set  the  boundary  conditions  at  infinity. 
Therefore,  the  Sommerfeld  radiation  conditions  have  to  be  replaced  by  some  other  boundary  conditions  set 
at  a  finite  location.  To  preserve  the  physics  of  the  model  that  involves  the  propagation  of  waves  toward 
infinity,  one  may  choose  to  set  the  appropriate  ABCs  at  the  external  artificial  boundary  that  truncates  our 
domain.^  We  emphasize  that  previously,  i.e.,  in  the  continuous  case,  we  have  introduced  and  used  the  ABCs, 
see  Section  2.2,  only  for  the  purpose  of  obtaining  compactly  supported  controls,  and  the  ABCs  applied  to 

^The  ABCs  that  would  guarantee  that  the  interior  solution  can  be  extended  beyond  the  artificial  boundary  so  that  the 
extension  solve  the  Helmholtz  equation  and  display  the  correct  far-field  behavior. 
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Fig.  3.1.  Schematic  geometry  of  the  domains,  the  stencil,  and  the  grid  boundary  7  =  7+  U  7  ;  Hollow  bullets  —  7+, 
filled  bullets  —  7“ . 


the  auxiliary  function-parameter  see  (2.13)  or  (2.17).  In  the  discrete  case,  the  ABCs  should  apply 

to  the  actual  solutions  ^  and  which  represent  the  acoustic  fields  on  the  grid.  For  the 

purpose  of  constructing  the  controls,  however,  we  will  never  need  to  implement  the  ABCs  for  the  acoustic 
solutions  on  the  grid  explicitly.  We  will  only  need  to  know  that  these  boundary  conditions  can  be  obtained 
(a  variety  of  different  approaches  can  be  found,  e.g.,  in  [17]),  and  that  they  will  guarantee  the  solvability  of 
the  difference  equations  involved. 

There  will,  of  course,  be  an  explicit  role  for  the  ABCs  in  the  discrete  framework  as  well.  We  will 
employ  these  boundary  conditions  in  the  same  capacity  as  we  have  used  the  original  continuous  ABCs 
in  (see  Section  2.2).  Namely,  when  obtaining  compactly  supported  controls,  the  ABCs  will  truncate  the 
corresponding  discrete  function-parameter  A  specific  approach  to  constructing  the  discrete  ABCs  that 
we  use  in  this  work  is  based  on  discretization  of  the  continuous  ABCs  of  Section  2.2,  and  we  discuss  it  in  the 
following  Section  3.3.  Altogether,  both  the  discrete  acoustic  fields  and  the  discrete-function  parameter  that  is 
used  for  constricting  the  control  sources,  see  formula  (3.3)  below,  are  going  to  satisfy  the  same  ABCs.  This  is 
similar  to  the  continuous  case,  when  both  the  solution  itself  (i.e.,  acoustic  field)  and  the  function-parameter 
w  satisfied  the  same  Sommerfeld  radiation  conditions  at  infinity. 

The  general  solution  for  the  discrete  control  sources  that  eliminate  the  unwanted  noise 
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on  N+  is  given  by  the  following  formula  (cf.  formula  (2.1)): 


mGM“ 


5 


(3.3) 


where  =  Wn\  n  G  N“,  is  a  special  auxiliary  grid  function-parameter  that  parameterizes  the  family 
of  controls  (3.3).  The  requirements  that  this  function  must  satisfy  are,  again,  rather  “loose,”  and  can 
be  considered  natural  discrete  counterparts  of  the  corresponding  requirements  of  the  continuous  function- 
parameter  ^^;(^c),  see  the  discussion  around  formula  (2.2)  in  Section  2.1.  Namely,  at  the  grid  boundary  7  the 
function  has  to  coincide  with  the  overall  acoustic  field  to  be  controlled: 


(h)  ^  .Ah) 

n  ji^ry  n  nG7’ 


(3.4) 


We  note  that  since,  e.g.,  for  the  second-order  discretizations  the  grid  boundary  7  contains  two  layers  of  nodes, 
7+  and  7“,  see  Figure  3.1,  then  specifying  the  corresponding  nodal  values  on  7  is  in  some  sense  equivalent 
to  specifying  the  function  and  its  normal  derivative  on  F  in  the  continuous  case,  see  (2.2).  Of  course,  this  is 
not  a  rigorous  statement  from  the  standpoint  of  approximation;  we  will  address  the  approximation-related 
issues  later  on.  We  also  note  that  when  creating  practical  designs,  the  boundary  data  shall  be 

interpreted  as  measurable  quantities  that  provide  input  for  the  control  system.  In  other  words,  we  can  think 
of  a  microphone  at  every  node  of  7;  these  microphones  measure  the  characteristics  of  the  actual  acoustic 
field  and  generate  the  input  signal 

The  other  requirement  of  besides  the  interface  boundary  conditions  (3.4),  has  already  been  men¬ 
tioned.  The  function  must  satisfy  the  appropriate  discrete  ABCs  at  a  finite  external  artificial  boundary. 
The  role  of  the  discrete  ABCs  is  the  same  as  that  of  the  continuous  ABCs  —  to  provide  a  replacement  for 
the  Sommerfeld  radiation  boundary  conditions.  This  is  done  in  the  same  approximate  sense  as  the  operator 
approximates  X,  see  Section  3.3.  Other  than  the  two  aforementioned  requirements,  the  function 
is  arbitrary  and  as  such,  parameterizes  a  substantial  variety  of  discrete  control  sources,  see  (3.3).  The  latter 
will  provide  the  search  space  for  optimization  in  Section  4. 

The  justification  for  formula  (3.3)  as  the  general  solution  for  the  discrete  control  sources  is  based  on 
the  theory  of  difference  potentials,  see  [13].  In  the  framework  of  this  theory  one  can  show  that  the  solution 
of  the  equation  subject  to  the  appropriate  ABCs,  where  is  defined  according 

to  (3.3)  and  (3.4),  will  be  equal  to  exactly  on  the  interior  sub-grid  N+ : 

In  other  words,  when  the  controls  g^^^  of  (3.3)  are  added  to  the  original  source  terms  of  the  governing  finite- 
difference  equation,  they  annihilate  the  unwanted  noise  on  the  domain  of  interest  in  the  discrete  sense,  i.e., 
on  the  grid.  The  aforementioned  solution  is  called  the  generalized  difference  potential  with  the 

density  —u^n^  defined  on  the  grid  boundary  7.  It  is  shown  in  the  theory  of  difference  potentials,  see  [13], 
that  the  potential  depends  only  on  its  density  and  not  on  the  values  of  the  function-parameter 

outside  7:  Consequently,  all  possible  controls  g^^^  obtained  according  to  (3.3)  with  different 

subject  only  to  (3.4)  and  the  corresponding  ABCs  (see  Section  3.3),  will  produce  identical  output 
on  N+  that  will  cancel  out  the  unwanted  noise  This  provides  room  for  optimization  of  the 

discrete  control  sources,  see  Section  4.  It  can  also  be  shown  that  every  discrete  control  source 
that  cancels  out  can  be  represented  in  the  form  (3.3)  with  some  function  n  G  N“, 

that  satisfies  (3.4)  and  the  external  boundary  conditions  (ABCs).  Similarly  to  the  continuous  case,  this 
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is  done  by  explicitly  constructing  the  appropriate  for  a  given  and  we  refer  the  reader 

to  [13,  Part  VIII]  and  [20,21]  for  detail. 

As  has  been  mentioned,  the  cancellation  of  noise  in  the  discrete  framework  is  obtained  on  the  grid  N+ . 
It  is  important  to  understand  in  what  sense  this  discrete  cancellation  models  the  continuous  cancellation 
described  in  Section  2.  This  is  basically  the  question  of  approximation  of  the  continuous  generalized  potentials 
by  the  discrete  ones.  To  that  effect,  the  theory  of  difference  potentials,  see  [13],  says  that  under  certain  natural 
conditions,  the  difference  potential  =  v^\  n  G  N+,  i.e.,  the  solution  to  with  g^^^  given 

by  (3.3),  approximates  the  continuous  Calderon’s  potential  v  =  v{x)^  x  e  ft,  see  (2.9),  i.e.,  the  solution 
to  ^  with  g  given  by  (2.1).  The  aforementioned  natural  conditions  include  first  the  consistency  and 

stability  of  the  finite-difference  scheme  for  the  Helmholtz  equation.  Consistency  and  stability  will  guarantee 
convergence  as  the  grid  size  vanishes.  In  addition,  the  discrete  boundary  data  of  (3.4)  have  to 

approximate  the  continuous  boundary  data  [u,  f^)  |p  of  (2.2)  in  the  following  sense.  Once  the  continuous 
function  u  and  its  first-order  normal  derivative  are  known  at  the  boundary  T,  normal  derivatives  of 
higher  orders  can  be  obtained  via  the  differential  equation  itself,  and  the  near-boundary  values  can 

be  calculated  using  Taylor’s  expansion;  the  order  of  accuracy  of  the  latter  calculation  with  respect  to  the  grid 
size  h  has  to  be  at  least  as  high  as  the  order  of  accuracy  of  the  interior  scheme.  In  this  case,  the  quality  of 
approximation,  i.e.,  the  rate  of  convergence  of  the  discrete  potential  to  the  continuous  one  with  respect  to  h, 
will  be  the  same  as  prescribed  by  the  finite-difference  scheme  itself.  For  the  second-order  central-difference 
schemes  discussed  in  Sections  3.3  and  4,  this  rate  is  0(/i^).  In  other  words,  when  designing  an  active  control 
system  following  the  finite-difference  approach,  one  can  expect  to  have  the  actual  noise  cancellation  in  the 
same  approximate  sense  as  the  solution  of  the  finite-difference  equation  approximates  the  corresponding 
solution  of  the  original  differential  equation.  Note,  in  any  particular  practical  setting  we  will  need  to  require 
sufficient  wave  resolution  on  the  grid,  i.e.,  the  waves  of  length  X  =  27r/k,  where  k  is  the  wavenumber  in  (1.1), 
will  have  to  be  well  resolved  by  the  specific  discretization. 


3.3.  Specific  Discretization  and  Discrete  Artificial  Boundary  Conditions.  As  has  been  men¬ 
tioned,  one  can  use  different  approaches  to  construct  discrete  ABCs  (see,  e.g.,  [17])  that  are  needed  to  obtain 
compactly  supported  controls.  The  most  straightforward  technique,  which  is  adopted  in  the  current  study, 
although  it  is  apparently  not  the  most  universal  one,  is  to  directly  approximate  the  continuous  boundary 
conditions  (2.13)  or  (2.17)  with  sufficient  order  of  accuracy.  For  that,  we  will  need  a  grid  that  would  be 
fitted  to  the  shape  of  the  external  artificial  boundary,  i.e.,  a  polar  or  spherical  grid.  An  example  of  the 
corresponding  grid  subsets  7+  and  7“  for  polar  coordinates  is  schematically  shown  in  Figure  3.2. 

In  all  numerical  experiments  that  follow  in  Section  4,  we  use  a  two-dimensional  setup.  Accordingly,  we 
introduce  a  polar  grid  that  has  J  cells  in  the  radial  direction  with  the  nodes  pj  =  j/Xp,  j  =  0, . . .  ,  J,  so 
that  y9o  =  0  and  pj  =  R]  and  L  cells  in  the  circumferential  direction  with  the  nodes  9s  =  s/X9,  s  =  0, . . .  ,  T, 
so  that  00  =  0  and  6l  =  For  simplicity,  it  is  convenient  to  assume  that  the  grid  sizes  Ap  =  RfJ  and 
AO  =  I'KjL  are  constant;  in  applications,  the  grid  in  the  radial  direction  may  be  stretched,  see  Section  4.2. 
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The  Helmholtz  equation  is  discretized  on  this  grid  with  the  second  oder  of  accuracy  by  central  differences: 


pj 


Ap 


*Xi=o. 


The  left-hand  side  of  (3.5)  is  a  particular  realization  of  the  operator  (3.1)  that  employs  the  five-node  stencil 
shown  in  Figure  3.2.  This  operator  will  be  used  in  Section  4  for  obtaining  optimal  discrete  control  sources. 


Fig.  3.2.  Schematic  geometry  of  the  domains,  the  stencil,  and  the  grid  boundary  7  =  7+U7  in  polar  coordinates:  Hollow 
bullets  —  7+;  filled  bullets  —  7“. 

To  construct  the  finite-difference  ABCs  at  y9  =  i?,  we  will  also  consider  a  semi-discrete  form  of  the 
homogeneous  equation  in  the  far  field: 

\  d  f  dWs  A  1  ^s+l  ~  Su’s  “h  Ws—i  .72  n  n  r  i  to 

z-rAp-rr]+-^ - - +  fcV  =  0,  s  =  0,...,L-l.  (3.6) 


p  dp  \  dp  J  p- 


Introducing  the  direct  and  inverse  discrete  Fourier  transforms,  /  =  — Z//2-hl,...  ,T/2,  s  =  0,...  ,Z/  —  1: 


L/2 

l^-L/2+l 
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we  reduce  (3.6)  to  the  following  system  of  ODEs  with  respect  to  wi  =  wi{p): 


1  d  ^ ^dwi^  af 


wi  +  k^wi  =  0, 


Oil  = 


4  .  2  lAe  ^  „ 

^sin  — ,  p>R, 


(3.8) 


p  dp  \  dp  J  p- 

where  again  I  =  —L/2  +  1, . . .  ,  I//2.  Equations  (3.8)  are  the  same  as  (2.14),  except  that  in  the  discrete  case 
the  range  for  I  is  finite,  and  f'  in  (2.14)  has  been  replaced  by  in  (3.8).  Therefore,  we  can  use  the  same 
boundary  conditions  (2.17)  for  I  =  —Lj2  +  1, . . .  ,  Lj2: 


d  . 

dp 


jH^hkR) 

P=R  H^ahkR) 


=  m{R) 


(3.9) 


only  with  the  Hankel  functions  of  order  I  replaced  by  the  Hankel  functions  of  the  order  a/ .  For  implementation 
in  the  foregoing  discrete  framework,  boundary  conditions  (3.9)  for  all  I  =  — Z//2  +  1,...  ,I//2  have  to  be 
approximated  with  the  second  order  of  accuracy,  which  can  be  easily  done  as  follows: 


^LJ  —  Wl  j-i 


Ap 


■A - ^ - =0, 


A 


_  dp 


(kR) 


H^ahkR) 


Finally,  relations  in  (3.10)  for  all  I  =  — Z//2  +  1, . . .  ,  I//2  can  be  rewritten  in  the  matrix  form: 

+  A  1  (  XT  “  A 


=  F  ^diag 


1 

Ap 


F  w.j-i  =  T  w.j- 


1: 


(3.10) 


(3.11) 


where  F  and  F~^  are  matrices  of  the  direct  and  inverse  discrete  Fourier  transforms  of  (3.7),  and  w.^j  and 
m,j_i  are  T-dimensional  vectors  of  components  and  respectively,  s  =  0, 1, . . .  ,  T  —  1. 

In  the  three-dimensional  case,  instead  of  the  discrete  Fourier  transforms  (3.7)  one  can  use  expansions 
with  respect  to  the  so-called  finite-difference  spherical  functions  (see  [13,  Part  IV,  Chapter  4])  that  form 
a  full  orthogonal  system  of  eigenvectors  for  the  spherical  part  of  the  discrete  Laplacian.  Other  than  that, 
the  construction  of  the  discrete  three-dimensional  ABCs  will  be  similar  to  the  foregoing  two-dimensional 
construction.  We  do  not  delineate  on  it  here  because  we  do  not  conduct  three-dimensional  computations 
in  this  paper.  Let  us  also  mention  that  for  small  I  the  difference  between  P  and  the  corresponding  see 
(3.8),  will  obviously  be  small  as  well.  Therefore,  for  smooth  functions  for  which  the  short-wave  part  of 
the  spectrum  (large  /’s)  is  insignificant,  one  may  not  even  have  to  replace  I  in  (2.17)  by  a/  in  (3.9).  We 
have  observed  this  type  of  behavior  in  the  previous  paper  [12],  in  which  we  studied  similar  questions  for  the 
Poisson  equation. 


3.4.  Types  of  Discrete  Control  Sources.  Similarly  to  the  continuous  case,  see  Section  2.3,  let  us 
now  identify  some  particular  types  of  the  discrete  control  sources.  First,  we  define  another  subset  of  the  grid 
M  (more  precisely,  of  M“): 

Wnt  =  {m  €  M-  |N„  n  7+  =  0}  . 


Basically,  is  the  interior  subset  of  M“ ,  such  that  when  the  center  of  the  stencil  sweeps  this  subset,  the 
stencil  itself  does  not  touch  7+,  see  Figures  3.1  and  3.2.  In  other  words,  we  can  say  that  is  a  subset  of 
M“  such  that 


U  Hn=^\7+- 

raeSCit 
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Having  defined  this  new  subset  introduce  the  auxiliary  function  =  w^\  n  G  N  ,  for  (3.3) 

as  follows: 


and 


w, 


ih)\ 


lnG7+ 


=  „wi 


lnG7+ ’ 


(3.12a) 


(h)  _ 

nG7“  nG7“  ’ 


=  0  on 


(3.12b) 


As  before,  we  also  assume  that  satisfies  the  appropriate  discrete  ABCs,  see,  e.g.,  (3.11).  Definition 
(3.12a)  means  that  on  the  interior  part  of  the  grid  boundary  7+  we  simply  set  equal  to  the  given 

Definition  (3.12b)  is  actually  a  discrete  exterior  boundary- value  problem  of 
the  Dirichlet  type.  Indeed,  everywhere  on  and  “outside”  the  exterior  part  of  the  grid  boundary  7“,  i.e.,  on 
N“  \7+,  the  grid  function  is  obtained  as  a  solution  of  the  homogeneous  equation  =  0  (enforced 

at  the  nodes  )  supplemented  by  the  boundary  data  on  7“ :  | |  ,  which  is  specified  for 

the  unknown  function  itself.  Note,  relation  (3.12a)  and  the  first  relation  (3.12b)  together  are  obviously 
equivalent  to  (3.4).  Therefore,  the  function  defined  via  (3.12a),  (3.12b)  falls  into  the  general  class  of 
used  for  obtaining  the  discrete  control  sources,  see  Section  3.2. 

Problem  (3.12b)  can  clearly  be  considered  a  finite-difference  counterpart  to  the  continuous  Dirichlet 
problem  (2.21).  Therefore,  its  is  natural  to  call  the  control  sources  =  ^^onTpoie  obtained  by  formulae 
(3.3),  (3.12a),  (3.12b)  the  discrete  surface  monopoles.  Indeed,  because  of  the  definition  of  given  by 
(3.12a)  and  (3.12b),  these  ^^onopoie  generally  speaking,  differ  from  zero  only  on  the  grid  set 
which  is  a  single  “curvilinear”  layer  of  nodes  of  grid  M  that  follows  the  geometry  of  T.  Accordingly,  the 
output  of  these  controls  can  be  called  the  discrete  single-layer  potential;  it  was  first  introduced  and  analyzed 
in  our  recent  paper  [18].  Let  us  emphasize  that  unlike  the  continuous  surface  monopoles  (2.24),  which 
belong  to  a  different  class  of  functions  rather  than  the  volumetric  sources  (2.1)  and  (2.2)  (singular  5-type 
distributions  vs.  regular  locally  integrable  functions),  the  foregoing  discrete  surface  monopoles  belong  to 
the  same  original  class  of  discrete  control  sources  (3.3)  and  (3.4).  They  can  be  considered  as  the  ultimate 
reduction  of  the  volumetric  discrete  controls  (3.3)  and  (3.4)  to  the  surface.  In  the  following  Section  4,  the 
discrete  surface  monopoles  will  play  a  fundamental  role  for  the  analysis  of  the  optimization  problems. 

Besides  the  discrete  surface  monopoles  and  the  corresponding  single-layer  potential,  one  can  also  define 
the  discrete  surface  dipoles  and,  accordingly,  the  double-layer  potential,  see  [18].  Grid  dipoles  are  introduced 
for  the  pairs  of  neighboring  nodes  so  that  the  nodes  in  the  pair  are  assigned  equal  in  magnitude  and  opposite  in 
sign  values.  The  control  sources  in  the  form  of  discrete  surface  dipoles  can  be  obtained  by  solving  a  special 
Neumann-type  discrete  exterior  boundary- value  problem  for  the  auxiliary  function  which  would  be 

analogous  to  the  continuous  problem  (2.22).  The  construction  of  surface  dipoles,  however,  is  somewhat  more 
elaborate  than  the  foregoing  definition  of  surface  monopoles.  And  because  in  this  paper  we  basically  focus 
on  the  monopole-type  sources  only,  we  are  not  going  to  further  delineate  here  on  the  issue  of  discrete  surface 
dipoles,  and  will  rather  refer  the  reader  to  our  paper  [18]  for  detail. 

4.  Optimization  of  Control  Sources.  Once  the  general  solution  for  controls  is  available,  in  either 
continuous  (2.1)  or  discrete  (3.3)  formulation,  the  next  step  is  to  decide  what  particular  element  of  this 
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large  family  of  functions  will  be  optimal  for  a  specific  setting.  There  is  a  multitude  of  possible  criteria  for 
optimality  that  one  can  use,  we  discuss  some  of  them  in  the  forthcoming  paper^.  We  should  also  emphasize 
that  in  many  practical  problems  the  cancellation  of  noise  is  only  approximate  and  as  such,  the  key  criterion 
for  optimization  (or  sometimes,  the  key  constraint)  is  the  quality  of  this  cancellation,  i.e.,  the  extent  of  noise 
reduction.  In  contradistinction  to  that,  in  this  paper  we  are  considering  ideal,  or  exact,  cancellation,  i.e., 
every  particular  control  field  from  either  the  continuous  (2.1)  or  discrete  (3.3)  family  completely  eliminates 
the  unwanted  noise  on  the  domain  of  interest.  Consequently,  the  criteria  for  optimality  of  the  controls  that 
we  can  employ  will  not  include  the  level  of  the  residual  noise  as  a  part  of  the  corresponding  function  of  merit, 
and  should  rather  depend  only  on  the  control  sources  themselves.  We  realize,  of  course,  that  at  a  later  stage 
of  the  work  we  will  also  need  to  look  into  the  issues  of  approximate,  rather  than  identical,  noise  cancellation, 
for  the  reason  of  further  reducing  the  costs.  In  this  case,  optimal  solutions  found  in  the  framework  of  the 
exact  cancellation  are  likely  to  provide  good  initial  guesses  for  subsequent  optimization  in  the  approximate 
framework.  Moreover,  it  will  probably  be  possible  to  use  some  results  from  the  approximation  theory  to 
deal  with  the  issues  of  approximate  noise  cancellation  once  we  have  solutions  for  the  exact  cancellation.  We 
expect  that  this  approach  will  be  much  faster  than  any  algorithm  of  combinatorial  type.  A  similar  reduction 
in  computational  complexity  of  optimization  was  outlined  in  our  earlier  work  [6]  on  the  optimal  distributed 
control  of  the  exterior  Stokes  flow. 

4.1.  Optimization  in  the  Sense  of  Li.  To  derive  a  meaningful  criterion  for  optimization  of  the 
control  sources,  let  us  first  briefly  discuss  the  physical  meaning  of  the  quantities  involved  in  the  formulation 
of  the  problem.  The  most  natural  way  to  interpret  the  held  variable  u  =  u{x)  (as  well  as  its  discrete 
counterpart  is  to  call  it  acoustic  pressure.  Indeed,  acoustic  pressure  is  the  quantity  which  is  directly 
measured  by  the  sensing  devices  (microphones),  and  as  such  can  be  immediately  supplied  as  the  required 
boundary  input  data  for  the  control  system,  see  formula  (2.2)  for  the  continuous  formulation  of  the  problem 
and  formula  (3.4)  for  the  discrete  formulation.  Then,  the  quantity  /  =  f{x)  that  represents  distributed 
monopole  sources  on  the  right-hand  side  of  the  original  Helmholtz  equation  (1.1),  as  well  as  its  discrete 
counterpart  ,  shall  be  interpreted  as  volume  acceleration  per  unit  volume  multiplied  by  the  ambient  fluid 
density  (see,  e.g.,  [8,9]).  Taking  into  account  that  we  are  only  considering  the  time-harmonic  case  so  that 
^(  • )  =  —iuj{  • ),  we  can  write 


f{x)  =  iugoqix),  (4.1) 

where  oo  is  the  frequency  of  the  original  temporal  oscillations,^  =  const  is  the  ambient  fluid  density,  and 
q{x)  is  the  volume  velocity  per  unit  volume.  Let  us  recall  [8,9],  that  the  actual  volume  velocity  is  defined  as 
the  integral  of  the  normal  component  of  the  fluid  particles’  velocity:  q  =  Jg  Unda,  evaluated  over  a  closed 
surface  S.  It  is  often  referred  to  as  the  acoustic  source  strength,  see  [8,9],  that  pertains  to  the  sources  located 
inside  S.  Accordingly,  the  volume  velocity  per  unit  volume  q  in  formula  (4.1)  is  given  hy  q  =  ^lim^y  Unda, 
where  V  is  the  volume  enclosed  by  S.  For  distributed  monopole  sources,  this  quantity  is  assumed  finite  and  is 
often  referred  to  as  the  acoustic  source  density  (i.e.,  acoustic  source  strength  per  unit  volume).  In  the  case  of 
an  isolated  point  monopole,  i.e.,  a  5-type  source,  which  can  be  represented  as  a  vanishingly  small  oscillating 

^  J.  Loncaric  and  S.  V.  Tsynkov,  Quadratic  optimization  in  the  problems  of  aetive  eontrol  of  sounds  in  progress,  2002. 

^We  remind  that  the  wavenumber  k  in  the  Helmholtz  equation  (1.1)  is  given  by  /c  =  a;/c,  where  c  is  the  speed  of  sound. 
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sphere  of  radius  e  with  surface  velocity  ?7e,  the  volume  velocity  can  be  introduced  as  ^  =  lim  When 

e — ^H“0 

the  strength  q  of  this  isolated  monopole  is  finite,  the  associated  source  density  is  formally  infinite  (which  is 
natural  to  expect  for  a  5-type  source).  In  the  case  of  a  continuous  distribution  of  sources,  the  relationship 
between  the  source  density  and  source  strength  is  standard  (like  between  the  mass  density  and  total  mass, 
electric  charge  density  and  total  charge,  etc.)  and  basically  says  that  the  integral  of  the  source  density  over 
a  given  region  is  equal  to  the  overall  source  strength  associated  with  this  region. 

Obviously,  the  physical  meaning  of  the  quantity  g  =  g{x)  of  (2.1)  that  describes  the  control  sources  is 
the  same  as  that  of  the  original  right-hand  side  /  =  /(cc)in(l.l).  If  the  field  variable  is  interpreted  as 
acoustic  pressure,  then  g{x)  shall  be  interpreted  as  the  volume  velocity  per  unit  volume  or  alternatively, 
the  acoustic  source  strength  per  unit  volume  (up  to  a  multiplicative  constant)  of  the  control  sources.  It  is 
important  to  mention  that  as  we  are  studying  the  time-harmonic  traveling  waves,  all  the  quantities  involved 
in  the  formulation  of  the  problem  are  complex- valued.  This  is  essential,  as  otherwise  it  would  have  been 
impossible  to  account  for  the  key  phenomenon  of  the  variation  of  phase  between  different  spatial  locations. 

Having  identified  the  physical  meaning  of  the  variables  involved  in  the  noise  control  model  that  we  have 
adopted,  we  would  argue  in  the  current  paper  for  selecting  the  optimal  control  sources  based  on  minimization 
of  their  overall  absolute  acoustic  source  strength.  Mathematically,  this  translates  into  the  minimization  of 
the  Li  norm  of  the  control  sources: 

llslli  =  j  \g{x)\dx — )-min,  (4.2) 

suppp 

where  the  search  space  for  minimization  in  (4.2)  includes  all  the  appropriate  auxiliary  functions  w{x),  by 
means  of  which  the  controls  g{x)  are  defined  (see  formulae  (2.1),  (2.2),  and  the  discussion  in  the  beginning 
of  Section  2.1).  The  advantage  of  using  this  criterion  for  optimization  is  that  it  has  a  clear  physical  inter¬ 
pretation.  Indeed,  the  total  absolute  volume  velocity,  i.e.,  the  integral  amplitude  of  surface  velocity  for  all 
the  sources,  is  to  be  driven  to  a  minimum,  which  can  be  a  characterization  of  the  corresponding  engineering 
devices.  For  comparison,  we  note  that  the  criterion  based  on  the  the  L2  norm: 


does  not  have  a  similar  clear  physical  interpretation,  although  as  indicated  below  (see  also  our  forthcoming 
paper,  footnote  4)  the  corresponding  numerical  optimization  problem  is  much  easier  to  solve.  Another 
advantage  of  using  the  Li  norm,  or  in  other  words,  the  overall  absolute  acoustic  source  strength,  as  the 
cost  function  for  optimization  (minimization),  is  that  it  characterizes  only  the  control  sources  themselves. 
This  is  a  convenient  distinction  compared,  e.g.,  to  the  power-based  criteria,  which  would  always  involve  the 
interaction  between  the  sources  and  the  field  they  operate  in.  This  interaction  is  often  referred  to  as  the 
“load”  on  the  sources  by  the  field,  see  [9],  and  may  lead  to  certain  types  of  degeneration  when  solving  the 
optimization  problem,  see  footnote  4. 

In  the  discrete  framework,  the  Li  minimization  problem  that  corresponds  to  (4.2)  is  formulated  as 
follows: 

||p('‘)||i=  (4.3) 


21 


where  Vm  accounts  for  the  volume  in  3D  or  area  in  2D  of  a  particular  grid  cell,  and  again,  the  search  space 
includes  all  the  appropriate  auxiliary  grid  functions  through  which  is  defined,  see  formula  (3.3). 
The  function  is  supposed  to  satisfy  boundary  conditions  (3.4)  on  the  interface,  and  the  selected  ABCs 
at  the  external  artificial  boundary;  for  the  two-dimensional  examples  analyzed  in  the  following  Section  4.2 
the  latter  will  be  boundary  conditions  (3.11). 

Let  us  now  return  to  the  definition  (3.3)  of  the  discrete  control  sources  and  adopt  the  polar  framework 

of  Section  3.3.  The  finite-difference  operator  can  obviously  be  interpreted  as  a  matrix  with  N  columns 
and  M  rows,  where  N  is  the  number  of  nodes  n  =  (s,  j)  of  the  grid  N“  such  that  the  corresponding  radial 
coordinate  pj  <  i?,  i.e.,  j  <  J,  and  M  is  the  number  of  nodes  m  =  (s,  j)  of  the  grid  M“  such  that  the 
corresponding  radial  coordinate  pj  <  i?,  i.e.,  j  <  J— 1.  Denote  by  w  the  vector  of  N  components 
such  that  n  eW~  and  j  <  J.  The  components  of  w  can  obviously  be  arranged  in  a  particular  way  so  that 
this  vector  can  then  be  decomposed  into  four  sub- vectors: 

w  =  [Wj,  Wo,  (4.4) 

where  Wj  contains  all  those  and  only  those  for  which  n  e  7,  w.^j  and  w.^j-i  correspond  to  the 
outermost  and  second  to  last  circles  of  the  polar  grid,  respectively,  as  in  formula  (3.11),  and  Wq  contains 
all  the  remaining  components  of  ^t7.  In  accordance  with  (4.4),  the  matrix  can  be  decomposed  into  four 
sub-matrices: 


lW  =  [A,B,C,D],  (4.5) 

where  the  number  of  rows  in  all  four  is  the  same  and  equal  to  M,  A  has  as  many  columns  as  there  are  nodes 
in  7  (we  denote  this  number  |7|,  C  and  D  each  has  L  columns  (see  Section  3.3),  and  the  number  of  columns 
in  B  is  obviously  N  —  jyj  —  2L. 

Using  representations  (4.4)  and  (4.5),  one  can  rewrite  the  optimization  problem  (4.3)  as  follows: 

\\V{Awry  -h  Bwo  +  Cw.^j-i  -h  Dw.^j)\\i  — min,  (4.6) 

where  the  norm  in  (4.6)  is  a  conventional  Li  norm  on  complex  M-dimensional  vectors,  and  U  is  an  M  x  M 
diagonal  matrix  with  the  entries  given  by  the  corresponding  cell  areas  Vm-  Next,  we  recall  that  the  search 
space  for  optimization  (4.3)  is  composed  of  all  the  appropriate  grid  functions  which  means  that  the 
vector  w  in  the  optimization  formulation  (4.6)  is,  in  fact,  subject  to  a  number  of  equality-type  constraints 
that  come  from  the  interface  conditions  (3.4)  and  ABCs  (3.11).  More  precisely,  the  first  sub-vector  in 
(4.4)  is  known  and  fixed  because  of  (3.4)  and  we  can  rewrite  (3.4)  as  ^t7^  =  tty,  where  tty  is  given.  The  last 
sub-vector  w.^j  in  (4.4)  is  a  function  of  according  to  (3.11).  Therefore,  we  can  conclude  that  only  tjt7o 

and  contain  free  variables  that  provide  the  search  space  for  optimization,  and  as  such  rewrite  (4.6)  as 

min  \\V{Bwo  -h  (C  -h  DT)w.^j-i  -h  Awj)\\i  =  min \\Ez  -  /||i,  (4.7) 

W0,W.^J-1  z 

where  E  =  V[B,  C  -h  DT]  is  an  M  x  (A^  —  |7|  -  L)  given  matrix,  z  =  [tt7o,  is  an  (A^  —  jyj  —  L)- 

dimensional  vector  of  unknowns,  and  f  —  —  VAw^  is  an  M-dimensional  known  vector  of  the  right-hand 
side.  Minimization  problem  (4.7)  is,  in  fact,  a  problem  of  finding  a  weak  solution  in  the  sense  of  Li  of  an 
overdetermined  complex  linear  system  Ez  =  f. 


22 


Let  us  first  note  that  the  most  conventional  weak  formulation  for  an  overdetermined  system  Ez  =  f 
would  be  that  in  the  sense  of  I/2,  rather  than  (4.7).  The  L2  minimization  problem:  \\Ez  — /II2  — >  min,  has 
proven  easy  to  solve  numerically  even  for  rather  complex  geometries.  It  does  not  require  the  Moore-Penrose 
type  arguments  and  can  be  conveniently  solved  by  a  standard  QR  algorithm;  we  report  the  corresponding 
results  in  our  forthcoming  paper  (see  footnote  4).  As  has  been  mentioned,  though,  this  formulation  lacks 
a  convincing  physical  interpretation  and  therefore,  hereafter  we  concentrate  on  solving  the  Li  optimization 
problem  (4.7). 

By  introducing  M  additional  real  variables  G  M,  i  =  1, . . .  ,  M,  one  can  reduce  problem  (4.7)  to  the 
following  optimization  problem  with  equality- type  constraints: 


mm 


in  ti 


-ti  =  0,  i  =  1, . . .  ,  M, 


which  is  equivalent  to  the  problem  with  inequality-type  constraints: 


(4.8) 


min^ti, 


i 


E' 


i  =  l, 


,M. 


3 

If  all  the  quantities  involved  in  the  formulation  (4.7)  were  real,  then  problem 
equivalent  to  the  linear  programming  problem  (see  [19,  Chapter  12,  §4]): 


(4.9) 


(4.9)  would,  in  turn,  be 


mm 


i 

^  ^  ^  ^ij^j  —  i  =  1,  .  .  .  ,  Af, 


(4.10) 


which  nowadays  can  be  solved  efficiently  even  for  large  dimensions.  However,  complex  entries  in  z,  and 
/,  see  (4.7),  are  essential  in  order  to  account  for  traveling  waves,  so  we  will  actually  need  to  solve  a  nonlinear 
problem  (4.9)  rather  than  a  linear  problem  (4.10). 

The  most  obvious  disadvantage  of  optimizing  in  the  sense  of  Li  is  that  the  foregoing  problem  (4.9) 
appears  very  difficult  to  solve  numerically.  Besides  being  nonlinear,  the  constraints  are  obviously  non¬ 
smooth.  Moreover,  strictly  speaking  those  constrains  are  not  convex  either.  Indeed,  for  every  i  =  1, . . .  ,  M, 
the  inequality  Xlj  ^ij^j  ~U  defines  a  cone  in  the  space  of  variables  5R(  Xlj  ^ij^j  and  eijZjy 

As  we  are  only  considering  the  upper  half  of  the  cone,  this  set  is  geometrically  convex.  However,  algebraically 

I  |2 

the  function 


—  ij  of  variables  Zj^  U,  obviously  cannot  be  convex.  And  the  algebraic  convexity  (i.e.. 


positive  semi-definiteness  of  the  Hessian)  is  exactly  what  distinguishes  between  the  convex  and  non-convex 
programming  problems,  with  the  latter  being  substantially  more  difficult  to  treat  in  a  numerical  setting, 
see  [19,  Chapter  24].  Of  course,  problem  (4.9)  can  be  reformulated  so  that  the  constraints  will  become  truly 
convex: 


mm 


—  ti<0,  i  =  1, . . .  ,  M. 


(4.11) 
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However,  in  the  formulation  (4.11)  the  most  “harmless”  cost  function  that  one  can  think  of,  i.e.,  the  linear 
function  has  been  replaced  by  the  function  that  has  singular  derivatives  at  the  optimum. 

Experimentally,  we  have  observed  that  this  presents  even  more  severe  problems  for  a  numerical  optimizer. 

Altogether,  the  combination  of  non-linearity,  non-smoothness,  and  only  “marginal”  convexity  (optimiza¬ 
tion  over  cones)  makes  problem  (4.9)  a  serious  challenge  even  for  the  most  sophisticated  state-of-the-art 
approaches  to  numerical  optimization  —  the  approaches  that  are  typically  based  on  interior  point  meth¬ 
ods  [10, 19].  The  difficulties  are  further  exacerbated  by  the  large  dimension  of  the  grid  that  the  problem  is 
formulated  on.  Even  for  the  aforementioned  state-of-the  art  methods  the  maximum  number  of  constraints 
that  they  can  handle  is  typically  on  the  order  of  hundreds.  And  in  problem  (4.9),  the  number  of  constraints 
is  the  same  as  the  number  of  grid  nodes  M.  As  such,  one  can  easily  encounter  an  orders  of  magnitude 
difference  between  the  number  of  constraints  that  the  numerical  optimizer  will  handle,  and  the  number  of 
grid  nodes  that  will  make  the  formulation  of  an  active  noise  control  problem  practically  interesting.  This  is 
especially  true  for  the  three-dimensional  problems. 

In  spite  of  all  difficulties,  we  have  still  managed  to  obtain  numerical  solutions  in  two  space  dimensions 
for  some  simple  test  cases.  All  numerical  experiments  that  we  have  conducted  indicate  a  very  consistent 
behavior  of  the  Li  optimal  solution  for  control  sources.  It  happens  to  be  the  discrete  layer  of  monopoles 
on  the  surface  ^monopfoie  described  in  Section  3.4.  Recall,  this  solution  is  obtained  by  applying  formula 
(3.3)  to  the  auxiliary  function  defined  by  (3.12a),  (3.12b).  In  the  following  Section  4.2,  we  report  the 
corresponding  computational  results,  and  in  the  subsequent  Section  4.3,  we  provide  a  general  proof  of  the 
global  Li  optimality  of  this  surface  monopole  solution  in  the  case  of  one  space  dimension. 

4.2.  Numerical  Solution  of  the  Li  Optimization  Problem.  For  our  numerical  simulations,  we 
have  considered  the  simplest  possible  two-dimensional  geometric  setup,  with  the  protected  domain  Q  in  the 
form  of  a  disk  of  radius  r  =  1  centered  at  the  origin.  The  external  artificial  boundary  was  a  circle  of  radius 
i?  >  r,  as  in  Section  3.3.  As  such,  the  resulting  discrete  control  sources  were  concentrated  within  the  annular 
region  r  <  p  <  R. 

In  contradistinction  to  Section  3.3,  we  have  used  here  a  polar  grid,  which  was  stretched  in  the  radial 
direction.  This  allowed  us  to  keep  the  cell  aspect  ratio  constant.  The  grid  is  first  built  in  the  coordinates 
(lny9,0),  it  has  equal  square  cells  x  ^  X  constructed  on  the  rectangle  [— ^,lni?]  x  [0,27r].  Then, 

the  conformal  mapping  maps  it  onto  a  polar  grid  with  uniform  angular  spacing:  6s  =  sA6,  where 

A6  =  ^  and  s  =  0, . . .  ,  T,  so  that  =  0  and  6l  =  Stt,  and  nonuniform  radial  spacing  pj  =  exp  •  j) , 
j  =  —1, 0, . . .  ,  J,  so  that  p-i  =  exp  (— ^) ,  po  =  1  =  r,  and  pj  =  R.  It  is  convenient  to  define  the  grid  sizes 
in  the  radial  direction  as  Apj  =  pj  —  pj^i  =  exp  •  j)  —  exp  •  (j  —  1)) ,  j  =  0, . . .  ,  J.  The  Helmholtz 
operator  can  be  easily  approximated  on  this  new  nonuniform  grid  with  the  second  order  of  accuracy  using  the 
same  five-node  stencil  as  shown  in  Figure  3.2.  This  involves  little  changes  compared  to  the  approximation 
(3.5),  which  works  for  uniform  grids,  and  we  refer  the  reader  to  our  paper  [12]  for  detail.  The  discrete  ABCs 
(3.10)  or  (3.11)  do  not  change,  except  that  Apj  needs  to  be  substituted  instead  of  Ap. 

As  we  are  building  our  control  sources  outside  of  the  protected  region  ft  =  {(p^6)  \  p  <  r  =  1},  i.e.,  on 
fti  =  we  do  not  need  to  be  concerned  with  the  structure  of  the  grid  inside  ft.  For  our  constructions, 

we  will  only  need  to  use  one  grid  circle  inside  ft.  This  will  be  the  innermost  circle  j  =  —1.  The  second  to 
innermost  circle  j  =  0  already  represents  the  interface  T  =  dft  =  {{p,6)  \  p  =  r  =  1}.  Adopting  the  definition 
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(3.2)  of  the  grid  subsets  introduced  in  Section  3.1,  we  obtain: 

=  =  M-  =  {(p,-,0,)|O<i< J-l}, 

=  {{pM\j  =  -1,0},  N-  =  {{Pj,9,)  I  -  1  <  j  <  J},  (4.12) 

7  =  li  = -1,0},  7+ =  {(Pj,6»s)  li  = -1},  7- =  {(pj,6is)|i  =  0}. 

For  all  definitions  in  (4.12),  we  assume  s  =  0, . . .  ,  L  —  1. 

In  the  computational  experiments,  we  have  used  grids  with  four  times  the  number  of  cells  in  the  circum¬ 
ferential  direction  compared  to  the  radial  direction.  Specific  grid  dimensions  were:  L  =  32  and  L  =  48,  and 
accordingly,  J  =  7  and  J  =  11.  (Note,  as  j  =  —1, 0, . . .  ,  J,  the  number  of  cells  in  the  radial  direction  is  J+1.) 
The  wavenumber  k  in  the  Helmholtz  equation  (1.1)  was  chosen  k  =  0.5.  The  excitation,  i.e.,  the  acoustic 
field  that  drives  the  control  system,  was  taken  in  the  analytic  form  of  a  shifted  fundamental  solution, 
see  formula  (2.3a),  as  if  it  were  generated  by  the  point  source  6{x  —  xi)^  where  Xi  =  (y9cos0,y9sin0)  =  (5,0). 
We  reemphasize  that  our  approach  does  not  require  an  explicit  knowledge  of  the  exterior  sources  of  noise. 
We  only  need  this  function  as  a  sample  field  to  be  used  as  given  data  in  formula  (3.4). 

We  have  also  considered  another  case:  T  =  48,  J  =  9;  for  this  case,  we  have  selected  k  =  0.9.  The 
excitation  was  produced  by  two  point  sources,  d{x  —  Xi)  -h  S{x  —  a^),  where  Xi  =  (5, 0)  and  X2  =  (1, 2). 

The  matrices  and  vectors  involved  in  the  formulation  of  the  optimization  problem  (4.7)  were  constructed 
in  accordance  with  the  chosen  geometric  setup.  Namely,  the  dimension  of  see  (4.5),  was  M  x  N  = 
(L-J)  X  (Z/-(J-h2)),  the  dimension  of  A,  which  corresponds  to  the  variables  on  7,  was  M  x2-L  =  (L-J)  x  2-Z/, 
the  dimension  of  B  was  M  x  {N  —  AL)  =  (T  •  J)  x  (T  •  (J  —  2)),  and  the  dimension  of  either  C  01  D  was 
MxL  =  {L’J)xL. 

We  have  tried  several  numerical  approaches  for  solving  the  corresponding  minimization  problems  (4.9), 
starting  with  the  algorithms  available  as  a  part  of  the  standard  Optimization  Toolbox  in  MATLAB.  However, 
our  best  numerical  results  were  obtained  with  the  software  package  SeDuMi  by  J.  F.  Sturm.®  This  is  a 
numerical  algorithm  for  optimization  over  cones  [15],  it  employs  the  ideas  of  interior-point  methods,  and 
the  self-dual  embedding  technique  of  [23],  see  also  [11].  The  algorithm  allows  for  complex-valued  entries, 
which  is  very  important  in  our  framework,  and  also  for  quasi-convex  quadratic  and  positive  semi-definite 
constraints.  Of  course,  all  the  cases  that  we  have  been  able  to  compute  using  SeDuMi  (see  above)  can 
still  be  treated  only  as  simple  model  examples  on  the  scale  of  potential  applications  for  noise  control  (see 
Section  4.1).  However,  the  optimal  solutions  that  we  have  obtained  all  demonstrate  a  very  coherent  behavior 
that  we  discuss  below.  On  Figures  4.1(a),  4.2(a),  and  4.3(a)  we  plot  magnitudes  of  the  Li  optimal  solutions 
computed  with  SeDuMi  [15].  Let  us  also  note  that  SeDuMi  is,  in  fact,  a  rather  general  procedure,  and 
one  may  expect  better  numerical  performance  from  more  focused  algorithms,  such  as  the  one  proposed  by 
Andersen,  et.  al.  in  [1].  In  the  future,  we  plan  on  trying  the  algorithm  of  [1]  for  solving  the  foregoing  Li 
minimization  problem. 

Apparently,  the  most  obvious  observation  that  one  can  make  by  looking  at  Figures  4.1(a),  4.2(a),  and 
4.3(a)  is  that  in  all  cases  the  optimal  solution  (i.e.,  the  Li  minimum)  is  concentrated  on  a  single  circumfer¬ 
ential  layer  of  grid  nodes.  This  is  the  second  to  innermost  circle  of  the  grid  N“ ,  see  (4.12),  i.e.,  the  grid  line 
j  =  0.  It  corresponds  to  the  outer  portion  of  the  grid  boundary  7“,  see  (4.12),  and  in  the  continuous  case  — 

®http : / / f ewcal . kub . nl/ sturm/ s  of tware/ sedumi . html 
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to  the  interface  T  itself,  F  =  \  p  =  r  =  1}.  In  other  words,  the  Li  optimal  solutions  for  control 

sources  that  we  have  computed  can  all  be  interpreted  as  layers  of  monopole  sources  on  the  perimeter  of  the 
protected  region  Q.  This  clearly  calls  for  comparing  these  optimal  solutions  with  the  densities  of  discrete 
single-layer  potentials  introduced  in  Section  3.4. 

To  conform  the  general  definitions  of  Section  3.4  to  the  specific  geometric  setup  analyzed  here,  we  need 
to  construct  on  N“  \7+  the  solution  of  the  discrete  exterior  Dirichlet  problem  (3.12b)  for  the  case  when 
7“  is  the  grid  circle  j  =  0  (in  this  case,  corresponds  to  j  >0).  Then,  the  operator  of  (3.5)  needs 
to  be  applied  to  the  overall  resulting  function  including  its  definition  (3.12a)  on  the  inner  part  7+  of 
the  grid  boundary:  .  Since  we  know  ahead  of  time  that  the  resulting  controls  will 

nG7+  nG7+ 

only  differ  from  zero  on  7“  (which  is  an  equivalent  of  in  this  case),  we  need  not  consider 

anywhere  beyond  j  =  1.  Consequently,  for  the  purpose  of  constructing  surface  monopoles,  we  may  simply 
set  J  =  1,  and  consider  on  three  grid  circles  only:  j  =  — 1,  j  =  0  =  J  —  1,  and  j  =  1  =  J.  In  so  doing, 
we  obviously  have  to  specify  the  ABCs  (3.11)  right  on  the  interface;  in  other  words,  the  input  for  the  ABCs 
will  be  on  7“,  .i.e.,  at  j  =  0,  and  the  output  —  on  the  outermost  circle  j  =  J  =  1.  Clearly,  specifying  the 
ABCs  on  7“  allows  us  to  reconstruct  for  j  =  J  directly  by  formula  (3.11),  i.e.,  without  actually  solving 
the  aforementioned  exterior  Dirichlet  problem.  And  once  we  know  for  j  =  —1,  0,  and  1,  we  can  easily 
obtain  the  discrete  surface  monopole  controls  on  7“,  i.e.,  for  j  =  0.  As  in  Section  3.4,  we  will  denote  these 
control  sources  SmoSoie- 


(a)  Li  optimal  solution  (b)  Surface  monopoles  ^^ono^oie 


Fig.  4.1.  Control  sources  for  L  —  32,  J  —  7,  k  —  0.5,  excitation  6{x  —  xi),  xi  —  (5,0). 
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(a)  Li  optimal  solution  (b)  Surface  monopoles  ^^ono^oie 

Fig.  4.2.  Control  sources  for  L  =  48,  J  =  11,  /c  =  0.5,  excitation  6{x  —  a;i),  xi  =  (5,0). 


(a)  Li  optimal  solution  (b)  Surface  monopoles  ^^ono^oie 

Fig.  4.3.  Control  sources  for  L  —  48,  J  —  9,  k  —  0.9,  excitation  S{x  —  xi)  +  5{x  —  X2)j  xi  —  (5,0),  X2  —  (I52). 

In  Figures  4.1(b),  4.2(b),  and  4.3(b),  we  plot  magnitudes  of  the  discrete  surface  controls  ^monopoie 
exact  same  cases,  for  which  we  have  explicitly  computed  the  Li  minimal  solutions  using  SeDuMi.  Comparing 
Figures  4.1(a),  4.2(a),  and  4.3(a)  with  respective  Figures  4.1(b),  4.2(b),  and  4.3(b)  at  a  glance,  we  conclude 
that  there  is  virtually  no  difference  between  them.  In  other  words,  the  Li  optimal  solutions  coincide  with 
the  surface  monopole  control  sources  ^mono^i^^oie-  corroborate  this  conclusion  accurately,  we  evaluate  the 
Li  norm  of  the  difference  on  the  grid  M“  between  each  Li  optimal  solution  and  the  corresponding  surface 
monopole  layer  ^^onopoie’  sissuming  that  ^^ono^poie  “  ^  everywhere  except  on  7“.  The  results  are  presented 
in  Table  4.1. 
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Table  4.1 

Comparison  of  the  computed  Li  optimal  solutions  with  surface  monopoles. 


Case 

l|5^^C!lelll.M- 

Ah)  _  (h,surfj  II 
"min  "monopole  1 1 1  ,M“ 

Relative  diff. 

Figure  4.1 

0.5764 

0.5761 

0.0067 

0.0117 

Figure  4.2 

0.5769 

0.5761 

0.0036 

0.0063 

Figure  4.3 

0.9760 

0.9750 

0.0083 

0.0085 

The  data  in  Table  4.1,  which  take  into  account  both  magnitude  and  phase,  do  corroborate  that  the 
respective  solutions  are  close  to  one  another.  Moreover,  by  comparing  the  second  and  third  rows  in  Table  4.1 
we  can  apparently  observe  the  phenomenon  of  grid  convergence.  Indeed,  the  case  of  Figure  4.2  is  computed 
on  the  grid  which  is  1.5  times  finer  in  each  direction  than  the  grid  of  Figure  4.1.  For  a  second-order  scheme, 
we  can  consequently  expect  a  drop  in  the  error  by  a  factor  of  ~  2.25,  which  we  indeed  see  in  Table  4.1  for  both 
the  absolute  and  relative  difference  between  the  Li  minimum  and  surface  monopoles  ^monopoie*  Even 
though  the  solutions  that  we  are  computing  are  obviously  not  smooth,  and  as  such,  the  grid  convergence  may 
be  difficult  to  justify  analytically,  the  foregoing  experimental  observation  certainly  makes  our  point  about 
the  coincidence  of  and  ^monopoie  more  convincing. 

Summarizing  the  foregoing  numerical  results,  we  can  see  that  in  all  the  cases  analyzed  the  minimum  for 
the  Li  norms  of  the  control  sources  g^^^  on  M“ ,  see  (4.3),  is  actually  given  by  the  Li  norm  of  ^monopoie* 

=  l|5monopolelll,M-- 

The  right-hand  side  of  the  previous  equality  can  be  recast  into  a  more  natural  form  by  noticing  that  surface 
controls  ^^onTpoie  defined  only  on  7“.  Then,  we  can  replace  the  Li  norm  over  the  two-dimensional  grid 
domain  M“  by  the  Li  norm  over  the  “one-dimensional”  grid  subset  7“.  This  will  bring  about  a  factor  of  Api 
because  obviously  the  cell  areas  Vm,  see  (4.3),  that  correspond  to  nodes  j  =  0  are  all  equal  and  proportional 
to  Api.  As  such,  we  obtain: 

““  II  i,M-  =  IkmonoSle  II 1,^-  '  (4-13) 

Equality  (4.13)  basically  conjectures  global  minimality  of  the  surface  monopole  solution  for  active  controls 
in  the  sense  of  Li. 

As  of  yet,  of  course,  we  can  only  claim  that  equality  (4.13)  holds  because  it  has  been  corroborated 
by  a  particular  collection  of  numerical  experiments  that  we  have  conducted.  However,  motivated  by  the 
consistency  of  our  experimental  observations  that  all  suggest  (4.13),  see  Figures  4.1,  4.2,  and  4.3,  and 
Table  4.1,  we  have  been  able  to  prove  a  general  result  on  the  global  Li  optimality  of  the  surface  monopole 
solution  for  controls  in  both  continuous  and  discrete  formulation  in  the  one-dimensional  case. 

4.3.  One-Dimensional  Proof  of  Global  Li  Optimality.  It  will  be  convenient  to  consider  simulta¬ 
neously  both  the  continuous  and  discrete  formulation  of  the  one-dimensional  noise  control  problem.  Let’s 
denote  the  independent  variable  x  G  M  and  introduce  the  one-dimensional  Helmholtz  equation  for  the  field 
variable  u  =  u{x)  (cf.  equation  (1.1)): 

Lu  =  ^  -h  k‘^u  =  f{x).  (4.14) 
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Then,  we  introduce  a  uniform  grid  =  n  •  /i,  n  =  0,  ±1,  ±2, . . . ,  of  variable  x,  and  approximate  equation 
(4.14)  with  the  second-order  central-difference  scheme 


u, 


(h) 

n+1 


-  2ui^^  +  ui% 


h? 


(4.15) 


Note,  we  are  using  here  the  same  subscript  “n”  for  both  the  discrete  unknown  function  and  discrete 
right-hand  side  because  for  the  particular  scheme  (4.15)  they  are  defined  on  the  same  grid.  We  will 
still  need  to  distinguish,  however,  between  the  grids  M  and  N  when  constructing  the  necessary  grid  subsets. 

Let  us  assume  that  our  protected  region  fl  corresponds  to  x  <  0  and  accordingly,  the  complementary 
region  fli  corresponds  to  x  >  0.  Then,  the  continuous  control  sources  will  be  given  by  (cf.  formula  (2.1)) 


a;>0 


(4.16) 


where  the  auxiliary  function  ^^;(x),  x  >  0,  is  supposed  to  satisfy  the  interface  conditions  (cf.  formula  (2.2)) 


^^;(0)  =  ^^(0), 


dw 


dx 


x—O 


du 

dx 


x—O 


(4.17) 


and  the  appropriate  ABC,  i.e.,  the  radiation  boundary  condition,  as  x  — >  -hoc.  The  quantities  ii(0)  and 
in  (4.17)  are  assumed  given. 

Next,  applying  the  definitions  of  Section  3.1  to  a  particular  stencil  given  by  (4.15),  we  will  have:  M+  = 
{m\m  =  n  =  —1,  —2, . . . },  M“  =  {m  |  m  =  n  =  0, 1, 2, . . . },  N+  =  {n\n  =  0,  —1,  —2, . . . },  N“  =  {n\n  = 
—1, 0, 1, 2, . . . },  and  7  =  {n  |  n  =  —1, 0}.  Accordingly,  the  discrete  one-dimensional  control  sources  will  be 
given  by  (cf.  formula  (3.3)) 


(h)  ^  ^n+l 
iJn 


- 


n>0 


(4.18) 


where  the  auxiliary  grid  function  is  supposed  to  satisfy  the  interface  conditions  on  7  =  {n\n  =  —1,0} 
(cf.  formula  (3.4)) 


w 


(h)  _  ,ih) 


-1 5 


W. 


(h)  _  ,ih) 


0 


^0 


(4.19) 


and  the  appropriate  ABC  at  infinity,  or  in  other  words,  for  large  n’s.  Again,  the  quantities  and  in 
(4.19)  are  considered  given. 

To  obtain  the  continuous  ABC,  we  assume  that  the  auxiliary  function  ^^;(x)  satisfies  the  homogeneous 
version  of  equation  (4.14):  Lw  =  0  for  x  >  X  >  0.  This  equation  has  two  linearly  independent  solutions: 
e-tkx  jg  ^  right-traveling  wave,  and  is  a  left  traveling  wave.  In  the  one-dimensional  framework,  we  obvi¬ 
ously  need  to  treat  the  right-traveling  wave  as  outgoing  for  the  artificial  outer  boundary  x  =  X.  Therefore, 
employing  the  same  mode  selection  principle  as  in  Section  2.2,  we  arrive  at  the  following  ABC  (cf.  formulae 
(2.13)  and  (2.17)) 

=  —ikw{X),  (4.20) 


dw 

dx 
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which  guarantees  that  only  one  of  the  two  aforementioned  linearly-independent  modes,  namely  will 

remain  in  the  composition  of  w{x)  for  x  >  X.  Note,  boundary  condition  (4.20)  can,  in  fact,  be  interpreted 
as  the  Sommerfeld  radiation  condition.  In  the  one-dimensional  case,  it  can  be  specified  at  a  finite  location, 
in  contradistinction  to  the  multidimensional  case  when  these  boundary  conditions  can  only  be  specified  at 
infinity,  see  formulae  (1.2a),  (1.2b).  Altogether,  the  continuous  auxiliary  function  w  =  w{x)  that  defines  the 
control  sources  g{x)  by  formula  (4.16),  is  specified  on  the  interval  [0,X]  and  satisfies  boundary  conditions 
(4.17)  and  (4.20). 

To  obtain  the  discrete  one-dimensional  ABC,  we  will  not  approximate  (4.20)  with  finite  differences,  as 
we  did  in  Section  3.3  for  the  polar  case,  when  it  was  basically  the  only  option.  We  will  rather  use  a  genuine 
finite-difference  approach,  which  has  shown  efficient  in  many  cases  (see  the  review  [17]),  and  which  was 
studied  in  our  recent  paper  [4]  for  a  more  complex  formulation  that  involves  a  high-order  approximation  to 
the  Helmholtz  equation.  Let’s  assume  that  the  auxiliary  grid  function  satisfies  the  homogeneous  version 
of  equation  (4.15):  |^  =  0  for  n  >  A”  >  0  (one  may  think  that  X  =  (A  —  1)  •  h).  This  homogeneous 

finite-difference  equation  has  two  linearly  independent  solutions:  q'^  and  q~'^^  where  q  and  q~^  are  roots  of 
the  corresponding  algebraic  characteristic  equation 

q^  -{2-k^h^)q  +  l  =  0.  (4.21) 


These  roots  are  given  by  the  formulae: 

^  =  1  —  “  ikh\j  1  —  q~^  =  1  —  +  ikh]J  1  —  (4.22) 

It  is  easy  to  see  from  (4.22)  that  for  small  h  the  discrete  wave  q'^  approximates  the  continuous  right-traveling 
wave  and  the  discrete  wave  q~'^  approximates  the  continuous  left-traveling  wave  Therefore, 

the  solution  q^  shall  be  interpreted  as  a  discrete  outgoing  wave,  and  q~'^  shall  be  interpreted  as  a  discrete 
incoming  wave,  for  the  external  artificial  boundary  n  =  N.  To  guarantee  the  radiation  of  waves,  we  need 
to  select  q'^  and  prohibit  q~'^^  or  in  other  words,  require  that  =  c  •  for  n  >  A,  where  c  =  const. 
Accordingly,  we  arrive  at  the  following  discrete  ABC: 


w 


(h) 

N 


(h) 


(4.23) 


which  guarantees  that  only  one  of  the  two  aforementioned  linearly  independent  solutions,  namely  q'^^  will 
remain  in  the  composition  of  for  n  >  N.  Altogether,  the  auxiliary  grid  function  that 

defines  the  control  sources  by  formula  (4.18)  is  specified  on  the  grid  subset  {n  |  n  =  —  1, 0, 1, . . .  ,  A}  and 
satisfies  boundary  conditions  (4.19)  and  (4.23). 

From  now  on,  we  will  be  considering  only  the  situation  with  no  interior  sources.  In  other  words,  the 
only  field  present  in  the  model  before  control  will  be  the  incoming  field  with  respect  to  the  protected  region 
=  {x  G  M|  X  <  0}.  In  the  continuous  case  it  can  be  expressed  as  u{x)  =  u~{x)  =  and  in  the  discrete 

case  —  as  =  u^n^~  =  Aq~'^^  where  A  =  const.  This  restriction,  in  fact,  presents  no  loss  of  generality. 
Indeed,  if  we  had  both  components,  u{x)  =  u~{x)  -h  ii+(x)  =  Ae'^^^  -h  and  chosen  w{x)  according 

to  (4.17)  and  (4.20),  then  we  could  have  replaced  this  w{x)  by  w{x)  =  w{x)  —  =  w{x)  —  u'^{x). 

The  function  w{x)  would  satisfy  the  new  interface  conditions  w{0)  =  t^“(0), 

of  (4.17),  and  the  same  original  ABC  (4.20).  Most  important,  the  control  sources  generated  by  this  new 
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auxiliary  function  according  to  (4.16)  will  be  the  exact  same  control  sources  as  those  generated  by  w{x): 
Lw  =  L[w-\-Be~'^^^]  =  Lw^  x  >0.  Similarly  in  the  discrete  case,  if  we  had  =  Aq~^-\-Bq^, 

then  the  auxiliary  functions  and  —  Bq^  =  would  generate  the  exact  same 

discrete  control  sources  according  to  (4.18).  Of  course,  the  foregoing  argument  is  in  complete  agreement 
with  the  general  discussion  of  Section  2  on  insensitivity  of  the  control  sources  to  the  interior  sound. 

In  the  continuous  one-dimensional  case,  the  interface  between  the  protected  region  =  {x  G  M|  x  <0} 
and  its  complement  fli  =  {x  G  M|x  >  0}  is  obviously  one  point  x  =  0.  To  construct  the  corresponding 
“surface”  monopole  controls,  we  consider  a  special  form  of  the  auxiliary  function  ^^;(x).  Namely,  if  the 
original  field  to  be  controlled  is  the  left-traveling  wave  that  propagates  into  ft,  u{x)  =  u~{x)  = 

X  <  0,  then  we  take  ^^;(x)  in  the  form  of  the  right-traveling  wave:  ^^;(x)  =  Ae~'^^^,  x  >  0.  Obviously,  this 
function  ^^;(x)  solves  the  homogeneous  equation  on  fti:  Lw  =  0,  x  >  0,  and  satisfies  the  ABC  (4.16).  It 
also  satisfies  the  Dirichlet  boundary  condition  at  the  interface  x  =  0:  ^^;(0)  =  ii(0).  Altogether,  we  see  that 
^^;(x)  selected  this  way  solves  the  one-dimensional  counterpart  of  the  exterior  Dirichlet  problem  (2.21)  that 
we  used  in  Section  2.3  to  obtain  the  control  sources  in  the  form  of  surface  monopoles.  According  to  the 
analysis  of  Section  2.3,  surface  monopoles  are  obtained  by  applying  the  operator  —X  to  the  function  v  of 
(2.19),  which  has  discontinuous  first  derivative  across  the  interface.  In  the  specific  one-dimensional  case  that 
we  are  studying  here,  this  function  is  given  by 


^^(x)  = 


Ae^^^ 


for  X  <  0 
for  X  >  0 


(4.24) 


Applying  the  operator  —X,  see  (4.14),  to  the  function  ^^(x)  of  (4.24)  in  the  sense  of  distributions,  see  [22], 
we  obtain  the  following  “surface”  (in  fact,  point)  monopole  control  source  (cf.  formula  (2.24)): 

5Spoie  =  2^*fc^(^)-  (4-25) 

To  obtain  the  discrete  “surface”  monopoles  in  the  one-dimensional  case,  we  need  to  consider  =  Aq~'^ 
for  n  <  0,  and  =  Aq^  for  n  >  0;  we  also  set  The  aforementioned  solves  the  discrete 

homogeneous  equation  =  0  forn  >  0,  satisfies  the  discrete  ABC  (4.23),  and  the  interface  conditions 

=  u^\  or  equivalently  (4.19).  In  other  words,  the  selected  solves  the  one-dimensional  version  of 
the  exterior  Dirichlet-type  problem  (3.12b),  (3.12a),  that  we  used  in  Section  3.4  to  derive  surface  monopole 
controls  in  the  discrete  framework.  Applying  the  operator  —L^^\  see  (4.15),  to  the  foregoing  function 
we  obtain  the  discrete  surface  control  source 


(/i,surf)  _ 

^monopole 


_^(2^+p) 

0 


for  n  =  0 
for  n  >  0 


(4.26) 


We  are  now  prepared  to  formulate  our  central  result  on  global  Li  optimality  of  surface  monopoles 
in  the  one-dimensional  discrete  framework.  For  any  function  that  satisfies  the  ABC  (4.23) 

and  the  interface  conditions  (4.19),  where  =  Aq~^,  ^  <  0,  the  Li  norm  of  the  corresponding  control 
sources  (4.18)  will  always  be  greater  or  equal  than  the  magnitude  of  the  surface  monopole  (4.26)  times 
the  grid  size  h:  11^ 


ih)\ 


1  > 


mm. 


{h,  surf) 
^monopole 


h. 


In  the  case  =  Aq'^  the  equality  is  achieved,  and  as  such: 


{h,  surf) 
^monopole 


h.  In  other  words,  the  following  theorem  holds. 
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Theorem  4.1.  Let  a  complex-valued  function  be  defined  on  the  grid  n  =  —1,0, 1, . . .  ,iV; 

where  N  >  0  can  be  arbitrary.  Let  where  A  ^  C  is  a  given  constant;  and  =  qWq^^  and 

=  Q'^n-i?  where  q  is  defined  by  formula  (4.22).  Then^ 


min 


N-l 

E 

n=0 


^’i+i  - 


(h)  ,  ^(h) 

«  '  ^n-1 


h? 


=  |A| 


(4.27) 


Proof.  Let  us  introduce  new  quantities  po,  Pi,...  ,  Pn-2  so  that  =  Po'Wq^\  =  piw^\  ..., 
=  Pn-2W^j!)^_2-  Then,  the  sum  on  the  left-hand  side  of  equation  (4.27)  can  be  recast  in  the  following 
form  (taking  into  account  that  =  qWq^^  and 


N-l 

E 

n—O 


'^n+l  -  +  wj-l 


/l2 

l^lbol 


PlP2  -  2pi  +  1  2 

-  +A; 


;i2 


+mipolbi|---bjv-3| 


PoPi-2po  +  l  ,  ,^2 


h? 


+  fc  Po 


+ 


P2P-i  -2p2  +  l  2 
-+k  P2 


h? 


+  ... 


PN-2Q  —  2pN-2  +  1  ,  7  2 

- h  k  Pn-2 


h^ 


Next,  we  introduce  new  notations:  po  =  q  zq^  pi  =  q  ^  zi^  ...,  pn-2  =  ^  +  zn-2:  where  q  is  defined  by 
(4.22),  and  all  the  quantities  are  generally  assumed  complex.  Using  these  new  notations,  we  can  rewrite  the 
generic  term  on  the  right-hand  side  of  the  previous  equality  as  follows: 


\MPo\\Pi  \  ■  ■  -  IPn 


Pn+lPn+2  2pjj-|-i  +  1  7  2 

- h  k  Pn+1 


h? 


\A\\q  +  Zo\\q  +  zi\ . .  .\q  +  z„ 


{q  +  z„+i){q  +  z„+2)  -  2{q  +  Zn+i)  +  1 


\A\\q  +  zo\\q  +  zi\ . .  .\q  +  z„ 


q^  2q  +  1  2  I  ^n+i^n+2  A  q{Zfi-^-i  +  Zfi-^2)  22^71-1-1  2 


h? 


+  k  q  + 


h? 


“h  k 


\^\\q  +  ^o\\q  Zi\ . .  .\q Zn 


mk  +  ^o|k  +  ^l|  ...  k  + 


^n+2{q  +  Zfi-\-i)  ^  ^  f  q  2  ^2 


h?  V  h‘^ 

^n-\-2{q  H“  Zfi-^i)  ^  Zfi-\-i 


+  k^ 


/i2 


In  the  last  chain  of  equalities,  expression  ^ — ^P^  +  k^q  turns  into  zero  by  virtue  of  the  characteristic  equation 
(4.21),  and  fi  =  q  —  2-\-  k^h^.  Using  the  definition  of  q  (4.22),  we  can  obtain: 


— 1  —  \k^k?  —  ikh\/l  —  \k‘^h?  +  k^h? 

—1  +  \k^k?  —  ikh\  1  —  \k‘^h? 

2  V  4 

2  V  4 

— 1  -h  f  ~  ^  ~  ^k^h^  -h  k^h^  —  ^k^h^  =  1. 
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Finally,  collecting  all  terms,  we  can  now  have: 
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E 
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—  +  2- 
7  0  I  ^ 


Z2(q  +  Zi)  ,  Z1 
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h?  h? 
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+  \A\\q  + zo\\q  + zi\...\q  + zisi-'i\ 


ZN-2 
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-fJ- 


> 


1^1 


7«-l  ,  7.2 


/l2 


+  A;" 


|A|  1^1  +  |A|  g|  1^1  - \A\\q  +  Zo\^-^  +  \A\\q  +  Zo\^-^\fi\ 


\A\\q  +  zo\\q  +  zi\^-^  +  \A\\q  +  zo\\q  +  zi\^-^\n\  -  . . .  +  \A\\q  +  zo\\q  +  zi\ . .  .\q  +  zn-sI  ^^^2  1^1  = 


1^1 


7«-l  ,  7,2 


+  k^ 


In  other  words,  we  have  obtained  the  inequality 
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E 
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W. 


(h) 

n+l 


-  2wi!'^  + 


h‘^ 


+  k^wi^^ 


(4.28) 


To  establish  the  result  of  the  Theorem,  i.e.,  formula  (4.27),  it  only  remains  to  show  that  there  will  be  a 
particular  =  w^\  for  which  inequality  (4.28)  transforms  into  the  equality.  Clearly,  the  equality  in 
formula  (4.28)  is  achieved  for  ^  >  0,  because  in  this  case  zq  =  zi  =  ...  =  zn-2  =  0.  This 

completes  the  proof.  □ 

Let  us  also  recall  that  if  we  multiply  the  sum  on  the  left-hand  side  of  either  formula  (4.27)  or  formula 

(4.28)  by  the  grid  size  /i,  we  obtain  the  discrete  Li  norm  of  the  control  sources  Therefore,  inequality 

(4.28)  transforms  into 


(/i,  surf) 
^monopole 


(4.29) 


and  consequently,  we  have,  in  effect,  demonstrated  the  global  Li  minimality  of  the  surface  control  sources: 


mm 


(/i,  surf) 
^monopole 


(4.30) 


Equality  (4.30)  is  a  one-dimensional  counterpart  of  (4.13),  but  unlike  the  experimentally  established  formula 
(4.13),  equality  (4.30)  has  been  proven  rigorously. 

The  foregoing  proof  of  Theorem  4.1  also  reveals  the  mechanism  of  discrepancy  between  the  optimal 
control  ^monopoie  Other  suboptimal  controls.  Namely,  every  time  the  auxiliary  function  “departs” 

from  the  pure  right-traveling  wave  which  is  equivalent  to  having  ^  0,  we  may  pick  up  additional 
value  of  11^  in  case  the  actual  estimate  based  on  the  triangle  inequality  that  we  use: 
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1-^71+21 
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happens  to  be  “strictly  greater”  rather  than  “greater  or  equal”  for  this  given  term.  It  is  also  interesting  to 
look  into  the  role  of  the  ABC  (4.23).  This  boundary  condition  “swallows”  the  last  term  in  the  sum  so  that 
all  the  previous  terms  can  cancel  one  another  in  pairs,  and  only  the  first  term  ^^onopoie  remain. 

Next,  we  will  analyze  the  continuous  case.  Let’s  assume  that  w{x)  is  a  regular  smooth  function,  which 
is  defined  on  the  interval  [0,X]  and  satisfies  boundary  conditions  (4.17)  and  (4.20).  Let  us  also  assume  that 
the  gird  function  =  w^\  n  =  0, 1, . . .  ,  A/”  —  1,  is  the  trace  of  w{x)  on  the  aforementioned  uniform  grid 
with  size  h:  =  ^(^n)  =  ^(^  *  h).  Note,  we  need  to  require  sufficient  smoothness  of  w{x)  in  order 

to  guarantee  the  consistency  of  the  finite-difference  scheme:  =  Lw  -h  0{h?)  (see  equations  (4.14) 

and  (4.15)).  Let  us  additionally  define  sind  in  accordance  with  the  boundary 

conditions  (4.19)  and  (4.23),  respectively,  like  in  the  formulation  of  Theorem  4.1.  Then,  for  small  h  we  can 
disregard  the  quadratic  terms  in  the  definition  of  q,  see  (4.22),  and  have  for  the  right  endpoint: 
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(h) 
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—ikw 
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which  is  obviously  an  approximation  of  the  continuous  ABC  (4.20)  with  the  accuracy  0{h).  Similarly,  for 
the  left  endpoint  we  obtain: 
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W  _  ,,Xh) 


'0 


—  w: 
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^  =ikw^\ 


which  is  an  0{h)  accurate  approximation  of  the  second  boundary  condition  (4.17)  under  the  assumption  that 
the  field  to  be  controlled  is  u{x)  =  u~{x)  =  x  <  0,  and  consequently,  ^  =  iku{G).  Altogether,  we 

have  constructed  a  grid  function  =  w^\  n  =  —  1, 0, . . .  ,  A”,  that  satisfies  the  conditions  of  Theorem  4.1 
and  also  approximates  on  the  grid  all  the  continuous  requirements  to  the  function  w  =  w{x). 

Let  us  now  again  multiply  both  sides  of  inequality  (4.28)  by  the  positive  quantity  /i,  and  consider 
independently  the  limit  on  its  right-hand  side,  and  the  limit  on  its  left-hand  side,  as  h  — >  -hO.  First,  we 
obtain: 


2^  +  P\h  =  \A\ 


h  — >  2\A\k,  as  h 


+0. 


Note,  the  limit  is  equal  to  the  magnitude  of  the  surface  monopole  in  the  continuous  formulation,  see  (4.25). 
Next,  on  the  left-hand  side  we  have: 


N-l 

E 

n—0 


{^)  I 


N-l 


n=0 


^^(Xn)  +  k'^w{Xn) 


h  +  Oik"^) 


X 


dx,  as  h  — >  -hO. 


Note,  the  limit  is  equal  to  the  Li  norm  ||^||i  of  the  continuous  control  sources  ^(x)  defined  by  formula  (4.16). 
As  inequality  (4.28)  holds  for  any  given  value  of  /i,  we  can  claim  that  it  will  also  hold  in  the  limit  h  — >  -hO. 
Therefore,  we  have  arrived  at  the  following 

Corollary  4.2.  Let  a  complex-valued  function  w  =  ^^;(x)  be  defined  on  [0,X].  Let  w{0)  =  A,  where 
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A  ^  C  is  a  given  constant;  and  w'{0)  =  ikw{0)  and  w'{X)  =  —ikw{X).  Then, 


+  k^w{x) 


dx  > 


(4.31) 


It  is  easy  to  see  that  the  requirements  of  w{x)  formulated  in  Corollary  4.2  are  equivalent  to  the  conditions 
that  guarantee  the  appropriateness  of  w{x)  for  constructing  the  control  sources  g{x)  using  (4.16),  see  formulae 
(4.17)  and  (4.20).  Therefore,  the  result  of  Corollary  4.2,  i.e.,  inequality  (4.31),  can  be  recast  as 

ll^lli  >  21^1^.  (4.32) 

On  the  right-hand  side  of  inequality  (4.32)  we  have  magnitude  of  the  “surface”  monopole  ^monopoie  defined 
by  formula  (4.25).  Note,  unlike  in  the  previously  considered  discrete  case,  when  the  minimal  solution 
^monopoL  (4.26)  was  an  element  of  the  same  class  of  control  sources  g^^^  defined  by  (4.18),  here  the 
minimum  ^monopoie(^)  defined  by  (4.25)  belongs  to  a  different  class  of  functions,  namely,  singular  (i.e.,  5- 
type)  distributions,  as  opposed  to  regular  (i.e.,  distributions.  In  other  words,  ^monopoie(^)  ^  ^1(^)5 

and  not  even  Z/i^^^^(M).  As  such,  we  cannot  introduce  the  Li  norm  of  ^monopoie(^)-  Therefore,  inequality 
(4.32)  formally  has  to  stay  the  way  it  is.  However,  symbolically  we  can,  of  course,  write 

f  |^^monopoleW|«^a;  =  /  \2Aik6{x)\ dx  =  2\A\k,  (4.33) 

Jr  Jr 

which  allows  us  to  “informally”  interpret  inequality  (4.32)  as  if  ^monopoie(^)  (4.25)  provided  a  lower  bound 
in  Li  for  all  the  control  sources  g{x)  defined  by  (4.16).  Let  us  also  note  that  “integration”  with  respect  to  x 
in  (4.33)  that  “removes”  the  5-function  itself  and  leaves  only  its  magnitude  2fc|A|,  is  a  continuous  analogue 
of  multiplication  by  h  on  the  right-hand  side  of  formulae  (4.29)  or  (4.30).  Therefore,  we  conclude  that  the 
continuous  inequality  (4.32)  is  a  direct  counterpart  of  the  discrete  inequality  (4.29). 

Even  though  ^monopoie(^)  ^  Ti^^^^(M),  we  will  still  show  that  there  are  regular  control  sources  g{x)  G 
^(ioc)(j^)  arbitrarily  close  to  ^monopoie(^)  ^he  weak  sense.  More  precisely,  we  will  construct  a 

sequence  of  regular  auxiliary  functions  w^{x),  such  that  the  corresponding  ge{x)  obtained  according  to  (4.16), 
will  converge  to  ^monopoie(^)  (4.25)  in  the  sense  of  distributions.  This  will  allow  us  to  claim  that  although 
the  minimal  solution  ^monopoie(^)  singular,  it  is,  in  fact,  “on  the  borderline”  of  the  class  of  regular  solutions. 
In  other  words,  it  is  a  limiting  point,  in  the  sense  of  weak  convergence,  of  the  space  of  all  g{x)  defined  by 
formula  (4.25).  In  additions,  we  will  also  show  that  the  Li  norms  ||^e||  converge  to  the  magnitude  2k\A\  of 
the  “surface”  monopole  ^monopoie(^)  (4.25).  This  will  allow  us  to  formulate  in  the  continuous  case  the 
result  similar  to  the  minimality  (4.30)  but  in  the  sense  of  “infimum”  rather  than  “minimum.” 

Consider  a  regular  function  =  w^{x)  that  is  defined  on  [0,X]  and  satisfies  boundary  conditions  (4.17) 
and  (4.20)  with  u{x)  =  Ae'^^^  for  x  <  0.  In  addition,  let  us  assume  that  not  only  for  x  >  X,  but  also 
in  between  some  (small)  e  <  X  and  X,  the  function  We{x)  already  coincides  with  a  right-traveling  wave: 
We{x)  =  e  <  X  <  X.  In  other  words,  we  require  that  ^^;e(0)  =  A,  ^^;'(0)  =  ikA,  and  w[{e)  = 

—ikw^{e).  For  the  purpose  of  obtaining  the  aforementioned  convergent  sequence,  we  will  subsequently  let 
e  — >  -hO.  The  control  sources  ^e(^)  are  defined  according  to  formula  (4.16):  ^e(^)  =  —w”{x)  —  k^We{x), 
0  <  X  <  e,  and  for  x  >  e  we  have  ^e(^)  =0.  If  ip  =  ip(x)  is  a  test  function  on  M,  i.e.,  a  compactly  supported 
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infinitely  smooth  function,  see  [22],  then  the  corresponding  functional,  i.e.,  the  distribution  itself,  can  be 
represented  as  follows: 


= 

Jr  J 

0 

e 

{e)ip'{e)  —  w{0)ip{0)  —  [ip{e)w'{e)  —  </p(0)^^;'(0)]  +  J  [—w{x)ip"  (x)  +  k^w{x)ip{x)]dx. 


As  w{x)  G  and  ip{x)  is  a  test  function,  the  integral  on  the  right-hand  side  of  the  previous  equality 

vanishes  as  e  — >  -hO.  Let  us  now  additionally  assume  that  the  functions  =  We{x)  are  constructed  so 
that  We{e)  — >  ^e(O)  when  e  — >  -hO.  In  other  words,  we  assume  continuity  at  x  =  0.  Then,  because  of  the 
boundary  condition  (4.20)  at  x  =  e,  we  have  w[{e)  — >  —ikA  as  e  — >  -hO.  Altogether,  we  obtain 


{9e:^)  — ^  2ikAip{0), 


as 


+0. 


(4.34) 


The  limit  (4.34)  implies  that  in  the  sense  of  distributions 

g,(x)  ‘?monopoie(»)  =  ^ikAd(x),  as  e  — ^  +0.  (4.35) 


Let  now  specify  a  particular  form  of  ^^;e(x): 

_ 'ikA. 

We{x)  =  - X^  -h  ikAx  -h  A,  xG[0,e].  (4.36) 

For  this  function,  we  have  ^^;e(0)  =  We{t)  =  A,  ^^;'(0)  =  ifcA,  and  u’'(e)  =  —ikA^  and  consequently,  ^^;e(x)  of 
(4.36)  meets  all  the  previous  conditions.  For  the  Li  norm  of  the  corresponding  control  ^e(x),  we  obtain: 

e  e 

llPelll  =  J  \w”ix)+k'^We{x)\dx  =  J  |A| 

0  0 

As  we  always  have  0  <  x  <  e,  the  dominant  term  in  the  last  integral  for  small  e  is  ^ ,  and  therefore: 


7  2  /  k^  2  7  9 

k^  -h  ( - x^  -h  k^x - 

e  e 


dx. 


11^. 


elll 


2fc|A|, 


as 


+0. 


(4.37) 


Putting  together  the  result  of  Corollary  4.2,  as  well  as  limits  (4.35)  and  (4.37),  we  arrive  at  the  following 
Theorem  4.3.  Let  a  complex-valued  function  w  =  ^^;(x)  be  defined  on  [0,X].  Let  w{0)  =  A,  where 
A  e  C  is  a  given  constant;  and  w'{0)  =  ikw{0)  and  w'{X)  =  —ikw{X).  Then,  in  the  class  of  regular 
control  sources  g{x)  defined  by  formula  (4A6)  for  all  such  w{x),  one  can  identify  a  sequence  ^e(x)  that  would 
converge  in  the  sense  of  distributions  to  the  point  monopole  ^monopoie(^)  formula  (4.25): 

9 Ax)  -5-  gtonopoiAx)  =  2ikA6(x),  as  e  — >  +0. 

Besides,  the  magnitude  of  the  point  monopole  ^monopoie(^)  (4-^^)  provides  the  greatest  lower  bound  for  Li 
norms  of  all  the  control  sources  g{x)  of  (4-16): 


inf  ll^lli  =  2fc|A|.  (4.38) 

w{x) 

Clearly,  estimate  (4.38)  can  be  interpreted  as  global  Li  minimality  of  the  “surface”  control  (4.25)  among  the 
continuous  one-dimensional  control  sources  (4.16).  It  is  a  continuous  version  of  the  previously  established 
discrete  result  (4.30). 
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5.  Discussion.  For  the  problem  of  active  control  of  sound,  we  have  systematically  described  time- 
harmonic  general  solutions  for  the  volume  and  surface  control  sources  in  the  continuous  and  discrete  formu¬ 
lation  of  the  problem.  These  control  sources  guarantee  the  identical  cancellation  on  unwanted  noise  on  a 
predetermined  region  of  interest.  We  have  also  proposed  a  criterion  for  optimization  of  the  resulting  control 
sources.  This  criterion  chooses  the  overall  absolute  acoustic  source  strength  as  the  cost  function  for  mini¬ 
mization,  and  as  such  admits  a  clear  physical  interpretation.  Mathematically,  it  translates  into  minimization 
of  complex- valued  functions  in  the  sense  of  Li ,  which  is  a  very  challenging  problem  from  the  standpoint  of 
numerical  implementation.  We  have  still  managed,  though,  to  compute  several  two-dimensional  numerical 
solutions  using  the  algorithm  SeDuMi.  All  these  solutions  demonstrate  a  coherent  behavior  —  the  mini¬ 
mum  is  achieved  on  the  surface  of  the  protected  region.  In  other  words,  the  minimum  is  delivered  by  the 
appropriate  surface  monopole  controls.  Therefore,  the  numerical  evidence  that  we  have  received  indicates 
that  surface  monopoles  may  provide  a  global  Li  minimum  for  the  control  sources  in  the  general  setting.  We 
have  been  able  to  rigorously  prove  this  result  in  the  one-dimensional  case  for  both  continuous  and  discrete 
formulation  of  the  problem.  Even  though  we  have  not  yet  been  able  to  prove  a  similar  result  for  a  general 
multi-dimensional  framework,  we  still  believe  that  it  is  true,  because  a  combination  of  the  two-dimensional 
numerical  evidence  and  a  one-dimensional  accurate  proof  cannot,  in  our  opinion,  be  a  mere  coincidence. 
Therefore,  we  put  it  forward  below  in  the  form  of  a  conjecture.  Let  us  remind,  that  according  to  (2.24)  the 
surface  monopole  controls  are  given  by 


5moiole(»)  =  - 


dw  du 
dn  dn 


SiT)  =  - 


dw  du 
dn 


dn 


5(r)  =  vix) 


■m, 


(5.1) 


where  w{x)  =  w{x)—u'^{x)  as  before,  and  w{x)  is  a  solution  to  the  exterior  Dirichlet  problem  (2.21).  Then, 
we  can  formulate  the  following 

Conjecture  5.1.  Let  a  complex-valued  function  w  =  w{x)  be  defined  on  and  let  it  be 

sufficiently  smooth  so  that  the  operator  L  of  (1.1)  can  be  applied  to  w{x)  on  its  entire  domain  in  the  classical 
sense,  and  the  result  Lw  be  locally  absolutely  integrable.  Let,  in  addition,  w{x)  satisfy  the  interface  conditions 
(2.2),  where  u  =  u{x)  is  a  given  field  to  be  controlled,  and  the  appropriate  Sommerfeld  radiation  boundary 
conditions  at  infinity,  (1.2a)  or  (1.2b).  Then,  the  greatest  lower  bound  for  the  Li  norms  of  all  the  control 
sources  g{x)  obtained  with  such  auxiliary  functions  w{x)  using  formula  (2.1),  is  given  by  the  Li  norm  on  T 
of  the  magnitude  of  surf  ace  monopoles  (5.1): 


inf 

w{x) 


/  \g{x)\dx  =  /  \v{x)\ds. 

JQi  Jr 


Alternatively,  we  can  rewrite  formula  (5.2)  as 


(5.2) 


inf  ||p(a:)||i,Qi  =  ||i/||i,r.  (5.3) 

W{X) 

Equality  (5.3)  is  a  multi-dimensional  generalization  of  (4.38).  Let  us  also  notice  that  equality  (4.13)  that 
was  obtained  on  the  basis  of  experimental  observations  in  two  space  dimensions,  can  be  considered  a  discrete 
two-dimensional  prototype  of  (5.3). 

In  the  formulation  of  Conjecture  5.1,  we  did  not  include  the  results  on  the  convergence  of  a  sequence 
of  volumetric  controls  to  the  surface  layer  iyS{T),  see  (5.1),  and  on  the  convergence  of  the  corresponding  Li 
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norms,  as  we  did  in  Theorem  4.3  for  the  one-dimensional  case.  We  believe,  though,  that  these  results  can 
be  easily  formulated  and  justified  in  the  multi-dimensional  framework  using  an  approach  similar  to  the  one 
that  we  have  used  in  the  one-dimensional  case.  The  key  missing  part,  however,  that  does  not  yet  allow  us 
to  transform  Conjecture  5.1  into  a  Theorem,  is  proving  that  surface  monopoles  provide  a  lower  bound  for 
the  volumetric  controls  in  the  sense  of  Ti,  whereas  showing  that  this  is  the  greatest  lower  bound  is  more 
straightforward.  The  analysis  of  this  problem  will  be  a  subject  of  our  future  research.  In  this  connection  we 
can  only  mention  that  at  least  in  the  two-dimensional  case  the  geometry  of  the  protected  region  fl  should 
not  be  a  limitation  when  constructing  a  general  proof.  If  one  can  prove  the  result  for  a  constant-width  linear 
strip  with  periodic  boundary  conditions  on  its  sides,  and  with  the  interface  T  being  a  segment  of  the  straight 
line  normal  to  the  sides  of  the  strip,  then  for  any  other  shape  the  same  result  can  likely  be  obtained  with 
the  help  of  a  conformal  mapping. 
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we  focus  on  optimizing  the  overall  absolute  acoustic  source  strength  of  the  control  sources.  Mathematically,  this 
amounts  to  the  minimization  of  multi-variable  complex-valued  functions  in  the  sense  of  Li  with  conical  constraints, 
which  are  only  “marginally”  convex.  The  corresponding  numerical  optimization  problem  appears  very  challenging 
even  for  the  most  sophisticated  state-of-the-art  methodologies,  and  even  when  the  dimension  of  the  grid  is  small, 
and  the  waves  are  long. 

Our  central  result  is  that  the  global  Li-optimal  solution  can,  in  fact,  be  obtained  without  solving  the  numeri¬ 
cal  optimization  problem.  This  solution  is  given  by  a  special  layer  of  monopole  sources  on  the  perimeter  of  the 
protected  region.  We  provide  a  rigorous  proof  of  the  global  Li  minimality  for  both  continuous  and  discrete  opti¬ 
mization  problems  in  the  one-dimensional  case.  We  also  provide  numerical  evidence  that  corroborates  our  result 
in  the  two-dimensional  case,  when  the  protected  domain  is  a  cylinder.  Even  though  we  cannot  fully  justify  it,  we 
believe  that  the  same  result  holds  in  the  general  case,  i.e.,  for  multi-dimensional  settings  and  domains  of  arbitrary 
shape.  We  formulate  it  as  a  conjecture  at  the  end  of  the  paper. 
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